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ABSTRACT 

The modern merger hypothesis offers a method of forming a new elliptical galaxy 
through merging two equal-mass, gas-rich disk galaxies fuelling a nuclear starburst 
followed by efficient quenching and dynamical stabilization. A key prediction of this 
scenario is a central concentration of young stars during the brief phase of morpholo¬ 
gical transformation from highly-disturbed remnant to new elliptical galaxy. To test 
this aspect of the merger hypothesis, we use integral field spectroscopy to track the 
stellar Balrner absorption and 4000A break strength indices as a function of galactic 
radius for 12 massive (M» ^ 10 10 Mq), nearby (z ^ 0.03), visually-selected plausible 
new ellipticals with blue-cloud optical colours and varying degrees of morphological 
peculiarities. We find that these index values and their radial dependence correlate 
with specific morphological features such that the most disturbed galaxies have the 
smallest 4000A break strengths and the largest Balrner absorption values. Overall, 
two-thirds of our sample are inconsistent with the predictions of the modern merger 
hypothesis. Of these eight, half exhibit signatures consistent with recent minor merger 
interactions. The other half have star formation histories similar to local, quiescent 
early-type galaxies. Of the remaining four galaxies, three have the strong morpho¬ 
logical disturbances and star-forming optical colours consistent with being remnants 
of recent, gas-rich major mergers, but exhibit a weak, central burst consistent with 
forming ~ 5% of their stars. The final galaxy possesses spectroscopic signatures of 
a strong, centrally-concentrated starburst and quiescent core optical colours indicat¬ 
ive of recent quenching (i.e., a post-starburst signature) as prescribed by the modern 
merger hypothesis. 
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1 INTRODUCTION 

IFhe stellar growth of high-mass (> 10 10 Mgh 2 ), quies¬ 
cent (non-star-forming) galaxies remains an important topic 
in galaxy evolution studies. Large and deep redshift sur¬ 
veys have conclusively shown that these red galaxies have 
grown significantly since at least redshift z ~ 2 — 3 in 
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terms of both their number density (e.g., 

Bell et al. 2004a: 

Brown et al. 20071: Martin et alJ 2007;: Brammer et al. 2009) 

and average sizes at a fixed stellar mass 

Khochfar & Silk 

2000; Trujillo et al.l 20061, 200?j; Lee et al. 

2013). The ob- 


served build up of this population and its conserved prepon¬ 
derance of early-type galaxy (ETG) morphologies at differ - 
ent epochs dBell et al] 2004a : lBlanton]l2006l : lBell et aI1l2012l ) 
requires a process (or sequence of processes) which both in¬ 
creases the number of galaxies with early-type morpholo¬ 
gies and quenches star formation. The currently accepted 
explanation for this phenomenon is the so-called blue-to- 
red migration ll Faber e t al.l 120071 : lEliche- Moral et, ~ahl l2010l : 
iGoncalves et all 20121 ) which posits that the transformation 
of star-forming (blue) disk galaxies into quiescent, red ETGs 
is responsible for the increasing population of red sequence 
galaxies over cosmic time. One possible, but well-accepted, 
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mechanism driving the bl ue-to-red migr ation is the modern 
merger hypothesis (iHopkins et al.l2008bl ). Under this model, 
a new massive quiescent elliptical galaxy is built by trans¬ 
forming two star-forming disk galaxies into one concentrated 
quiescent red ETG via a gas-rich major (i.e., nearly equal 
mass) merger. 

Yet, in detail, the red ETG population contains myriad 
morphological types in detail from early-type spiral disks, to 
lenticular (SO) ga laxies, to the well-known dichotomy of el¬ 
lipticals (e.g.. lKormendvll2009l) . making disentangling which 
migration processes are dominant a daunting task that likely 
depends on galaxy morphology, mass, environment, and red- 
shift. We have made progress in the local universe by fo¬ 
cusing on the red wa r d migrati on of pure-spheroid elliptical 
galaxies (iMcIntosh et al.ll2014l . hereafter, M14). This choice 
is motivated by the fact that the formation of new ellipticals 
is theoretically tied to a single me chanism that both trans¬ 
forms disks (e.g., major merging, iToomre fe Toomrelll972l : 
iBarnesI Il992l. and references therein) and is hierarchic ally 
motivated ( White fe ReesllT978l : iKauffmann et alll l 993l ) . 

Simulations of gas-rich major mergers provide a key 
prediction which allows us to qualitatively identify major 
merger remnants. During the merging process, gravitational 
torques deposit cold gas in the centre of the remnant and 
produce a centrally-concentrated burst of star formatio n 
jBarnes fc Hernauistlll99ll . 1 19961 : iMihos fc Hernauistlll996ll . 
According to the modern merger hypothesis, this is quickly 
followed by the ignition of an active galactic n ucleus (AGN) 
which ultimately quenches star fo rmation (ISpringel et al.l 
l2005i : ISchawinski et al 1l2007|, Ho n!) and forms a quiescent 
elliptical galaxy ( Hopkins et al] 20061 . l2008al l . This is sup¬ 
ported by the bluer colours of early-type galaxies host- 
in g an AGN, c ompared with their non-active counterparts 
dSanchez et al.ll20041 ). The merging process finally produces 
a relaxed remnant with the structure and ki nematics of ob ¬ 
served moderate- mass spheroidal galaxies (ICox et al.ll2006l : 
iNaab et al]|2006h . 

We, therefore, anticipate ellipticals which are plausible 
remnants of recent gas-rich mergers to have bluer colours 
as compared to their “normal” counterparts on the colour- 
mass relation due to recent, enhance d star formation . Such 
blue ellipticals are rare. For example, iFukueita et al.l (j2004j l 
found 2 of the 210 ellipticals in their sample of 1600 bright 
galaxies from the SDSS at z<0.12 to be blue ellipticals, 
M14 found 1602 blue ellipticals in their complete sample of 
63 ,454 ma ssive (M» > 10 1 o Mq1i - 2 ) galaxies at z<0.08, and 
iLee et ah! 120060 found 29 in their sample of 1949 galaxies 
with spectroscopic redshifts < 0.4 in the GOODS fields. Des¬ 
pite their rarity, these galaxies are important as they build 
the transition from star-forming disks in the blu e cloud to 
mass ive, quiescent ellipticals on the red sequence feell et al.1 
1200711 . The robust sample of blue elliptical galaxies seen in 
M14 provides a good starting point for testing methods to 
distinguish recent elliptical mass assembly. 

Over the last decade, spatially -resolved s p ectro - 
scopic surveys such as SAURON (IBacon et all l200ll l 
and ATLAS3D (ICappellari et all l201ll i have drastically 
enhanced our understanding of the content and structure 
of cluster an d field e llip ticals. L ik ewise f projects such 
as PINGS dR osales- Orteea et al.l l2010l l. TYPHOON 
( Sturch fe Madorel I 2 OI 2 I 1 . and the CALIFA survey 
( Sanchez et al.l ~ 20lj i~ have provided an unprecedented 


level of detail in spiral galaxies. Aside from serendipitous 
discoveries (e.g., iFukugita et al.l 1 20041) . little work has 
been brought forth from the population of blue elliptical 
galaxies. Further, previous studies of these galaxies have 
been limited to evalu a ting only the centr al regions (e.g., 
ISchawinski et all 120091 ; iToieiro et al] I2013T ) which provides 
a narrow perspective of the mass assembly history. In this 
work, we propose to analyse the 4000Abreak strength and 
Balmer absorption line indices out to ~ 3 half-light radii 
(covera ge previously only ach ieved for red elliptical galax¬ 
ies, e.g. iWeiimans et al.l I2009T) to examine the spatial star 
formation histories of this rare, but important population 
of elliptical galaxies in the local universe. 

Our paper is organized as follows. Section 2 details 
our sample selection, observation setups, and data reduc¬ 
tion pipeline. Section 3 outlines our model construction. 
In Section 4, we provide results from our comparison of 
the spatially-resolved spectra to a suite of star formation 
history (SFH) models. Section 5 contains the discussion 
of our results. Finally, in Section 6 we provide a sum¬ 
mary and conclusions. Throughout this paper we calcu¬ 
late comoving distances in the ACDM concordance cosmo¬ 
logy with = 0.3, Ha = 0.7, and assume a Hubble 
constant] of Ho = 70 km s^ 1 Mpc^ 1 . Magnitudes from 
the Sloan Digital Sky Survey (lYork et al.l l2000l . SDSS) 
are on the AB system such that jhab = m + Am 
where Am = (-0.036 +0.01 2, +0.010, +0.028, +0.040) for 
(u, g,r,i, z) (IYang et al.l 1200 7). 


2 DATA SAMPLE AND OBSERVATIONS 
2.1 Sample Selection 


To test the prediction of the modern merger hypothesis 
that gas-rich major mergers lead to the construction of 
massive ellipticals with temporarily blue colours, we ac¬ 
quire observationally-intensive IFU spectroscopy using long 
exposure times to obtain spatially-resolved spectra bey¬ 
ond the bright central bulge. We then use the H5 a and 
H7A Lick indices and the narrow-band 4000Abreak to 
map the recent star formation histories as a function of 
galactic radius for a collection of visually-classifiec] plaus¬ 
ible new ellipticals provided in M14. Using the New York 
University Value-Added Galaxy Catalogue (NYU-VAGC, 
iBlanton et al.l 20051) based o n the SDSS Data Release 4 
(DR4. lAdelman-McCarthv et al] 120061 ). M14 construct a 
complete sample of nearby (0.01 ^ z ^ 0.08), high-mass 
(M* ^ 10 10 Mq/i -2 ) galaxies. Using the rest frame (g — r) 
Pet rosian colo u rs, th e absolute r-band magnitudes, and 
the iBell et ahl |2003i ) M/L prescription to measure stel¬ 
lar mass, they apply the empirical colour cut 01 (g — r) ^ 
0.81 + 0.1 [log lo (M ga i i */M0h -2 ) — 10.0] to select blue-cloud 
galaxies (see Fig. [lj where the colors have been k-corrected 
to z = 0.1. They isolate early-type (spheroid-dominated) 


The stellar masses from M14 are calculated assuming Ho = 
lOO/i. 

The morphological classifications of M14 were determined by 
manual inspection of the images by people. Here and throughout 
the paper, we use the phrase “visually-classified” to mean classified 
“by eye.” 
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galaxies using th e common r- band concentration cut of 
R 90 /R 50 ^ 2.6 (IStrateva et al.1 [200 1:1. where Rgo (R 50 ) is 
the radius inside which 90% (50%) of the light is contained 
to produce a set of 8403 unusually blue ETGs. 

M14 then use visual classification with careful morpho¬ 
logical inspection to identify plausible new merger remnants 
from this set of blue ETGs. Visual inspection further al¬ 
lows them to remove contamination by bulge-dominated 
spirals (S) and inclined disks (iD), plus galaxies with un¬ 
certain morphology (U), and galaxies affected by artefacts. 
Their final catalogue consists of 1602 massive, blue elliptical 
(E,85%), peculiar elliptical (pE,8%), and spheroidal post¬ 
merger (SPM,7%) galaxies with high classifier-to-classifier 
agreement, defined as a minimum of three out of four classi¬ 
fiers in accord. The details of these identifications and their 
robustness are described in detail in M14. Briefly, the E, 
pE and SPM types represent three morphological bins that 
plausibly span a qualitative time sequence since merging (as¬ 
suming that all blue elliptical galaxies are derived from gas- 
rich major merging). At one extreme (E), galaxies appear 
to be dynamically relaxed ellipticals with little or no evid¬ 
ence of recent tidal activity, contrasted by those at the other 
extreme which appear to be freshly coalesced with very dis¬ 
turbed morphologies (SPM). In between the two extremes 
(pE), we find galaxies which appear to be relaxed elliptic¬ 
als with either a modest external morphological peculiarity 
(e.g., a loop or fan feature) or an internal structure like a 
dust lane. We parametrize the level of visual tidal features 
by defining a plausible post-merger type, T ppm , for each blue 
ETG based on the sum of its four classifications weighted as 
follows: 

4 

T ppm — ^ ) W c ,i (1) 

1=1 

where w c ,t is the i th classification such that w c = -1 (S/iD), 
-0.5 (U),'l (E), 2 (pE), and 3 (SPM); e.g., T ppm = 8 for 
four-way pE agreement. We use negative weights to dis¬ 
tinguish definite disks and uncertain classifications from 
(positive w c ) plausible post-merger systems. We investig¬ 
ate a subset of the 8403 blue ETGs from M14 using a re¬ 
laxed classifier agreement such that T ppm ^ 2.0 to produce a 
sample of 1915 E (71%), pE (6%), SPM (6%), plus 17% with 
combinations of E+pE+SPM classifications (i.e., classifiers 
agree that these objects are not S or iD despite disagree¬ 
ment over specific plausible post-merger type). For example, 
nyul00917 appears to be a relatively smooth elliptical with 
a large loop structure (see Figure [2]). The classifiers in M14 
were evenly split between the pE and SPM classification: 
falling below the 75% classifier agreement requirement. Us¬ 
ing our relaxed classifier allows this galaxy to be included in 
our sample. 

To maximize spatial coverage of the spectra, we require 
each galaxy to fill the FOV of its respective spectrograph 
(each spectrograph used in this work is~ 70” on a side) in 
one of two ways: (i) twice the galaxy’s R90 ~ 70” or (ii) the 
galaxy has loops, arms, or other structures which extend 
beyond its main body to fill the FOV. These requirements 
reduce the range of redshifts for our targets to 2 < 0.03 and 
provide a sample of 111 E (57%), pE (11%), SPM (10%), 
and other (T ppm ^ 2,23%) galaxies. Based on observational 
availability, we select a set of 12 (3 E, 3 pE, 2 SPM, and 4 
other) which nearly equally sample the morphology space 



log(M,/M 0 h 2 ) 


Figure 1. Colour-mass diagram for the 63454 nearby (0.01 ^ 2 Si 
0.08), high-mass ( M * 0 10 10 Mq/i — 2 ) SDSS DR4 galaxies from 
M14 with our sample of 12 blue ellipticals overlaid. The red line 
shows the M14 colour cut, symbols represent visual classifications 
from M14. Refer to Table [l] for letter labels. 

while also selecting interesting candidates for unique star 
formation histories based on a range of morphological details 
such as circumnuclear blue rings or large dust features. In 
Table [T] we list physical characteristics and morphological 
classifications for our sample. 

In Fig. |T] we show our sample in the colour-mass plane 
using the SDSS ( g-r ) colours k-corrected to z = 0.1 in the 
AB magnitude system and the log of the stellar mass in units 
of Mq1i™ 2 . With the exception of the two low-mass ellipt¬ 
icals with smooth morphologies, we find that the low-mass 
galaxies tend to have the strongest morphological peculi¬ 
arities and blue colours consistent with star-forming galax¬ 
ies. The more massive galaxies tend to exhibit green val¬ 
ley colours and smoother morphologies. We note that all 
of the peculiar ellipticals have masses between 10 10 ’ 5 and 
10 11 M 1 T 1 h~ 2 , above which extremely few blue galaxies lie 
(e.g., iBaldrv et al.ll2004l ). 

2.2 Observations 

Our final sample of 12 plausible merger remnants were 
observed over three separate runs between January 2010 
and April 2011. We require moderate resolution (R~1000) 
spatially-resolved spectra with wavelength coverage span¬ 
ning the 4000A break and Balmer Ha and H 7 absorption 
indices to obtain SFHs as a function of galactic radius to 
differentiate between remnants of major mergers from those 
of other mass assembly mechanisms. The optimal choice for 
large-field, spatially-resolved spectroscopy is an integral field 
unit (IFU) which additionally provides concurrent sky flux 
measurement to facilitate removal of contaminating sky lines 
during reduction. 

The first observing run was performed between Janu¬ 
ary 18 th and January 20 th , 2010 at the 3.5-meter Calar 
Alto observatory in Spain using the Po tsda m MultiAper- 
ture Spectrophotometer (PMA S, lRoth et al.ll2005h and the 
PPAK IFU llKelz et al.l |200 gIi . We used the V1200 grat¬ 
ing with first order and forward blaze angle centred on 
4263A, providing a spectral resolution FWHM of 2.7A, and 
an observed wavelength range of 3877A to 4633A. The 
spectra were binned 2x2, giving an instrumental resolu¬ 
tion of 0.75A/pixel. The 2.7” fibres provide a spatial res¬ 
olution of 1.63 kpc at z=0.03. Three thirty-minute expos- 
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Table 1. Sample 


Name 

(i) 

^published 

(2) 

R50 

(3) 

C r 

(4) 

g 

(5) 

M r 

(6) 

M* 

(7) 

Classification! 

(8) 

Agreement 

(9) 

Tppm 

(10) 

ID 

(11) 

nyul00917 

0.027 

2.9 

3.6 

15.10 

-21.11 

10.24 

pE/SPM 

50 

10 

.1 

nyu22318 

0.030 

5.9 

3.0 

14.79 

-21.72 

10.79 

pE 

75 

7 

G 

nyu535845 

0.014 

6.2 

3.0 

13.99 

-20.81 

10.24 

E 

75 

5 

D 

nyu541044 

0.020(“) 

6.0 

3.5 

13.55 

-22.24 

10.61 

SPM 

75 

11 

K 

nyu598180 

0.018 

4.4 

2.9 

14.58 

-20.75 

10.19 

E 

75 

5 

C 

nyu619567 

0.029 

5.2 

3.1 

14.16 

-22.31 

10.85 

E/pE 

50 

6 

F 

nyu674881 

0.022 

8.3 

3.3 

13.79 

-22.10 

10.80 

pE/U 

50 

3 

A 

nyu719486 

0.028 

2.4 

3.3 

15.08 

-20.85 

10.06 

pE/SPM 

50 

10 

I 

nyu742039 

0.025 

5.2 

3.1 

14.16 

-21.90 

10.70 

E 

100 

4 

B 

nyu835691 

0.026 (6 > 

5.3 

3.4 

14.42 

-21.41 

10.40 

SPM 

75 

11 

L 

nyu916578 

0.029 

5.1 

3.4 

14.85 

-21.34 

10.51 

pE 

75 

5.5 

E 

nyu916757 

0.028 

6.3 

3.1 

14.19 

-22.11 

10.77 

pE 

75 

7 

H 


Note; For each galaxy in our sample, we list (1) the ID provided in the NYU-VAGC from the SDSS DR4, (2) the spectroscopic redshift 
published in SDSS or NED (objects with a superscript are from NED), (3) the SDSS r-band Petrosian half-light radius in arcseconds, 
(4) the concentration index given by the ratio of the SDSS r-band Petrosian R90 (the radius inside which ninety percent of the light is 
contained) and the half-light radius (Rgo/Rso)? (5) the SDSS g-band apparent magnitude, (6) the extinction- and k-corrected (to z = 0 
absolute SDSS r-band magnitude, (7) the stellar mass in units of log (Mq h ~ 2 ), (8) morphological classification (see text for details), 
(9) classifier agreement percentage, (10) the plausible post-merger type parameter (Equation [TJ, and (11) the ID used to refer to each 

galaxy in figures ordered by increasing Tp pm . 
f All cl assifications are given where (9) i s less tha n 7 5 per- cent. 

(a) iRothberg Sz Joseohl l|2006ah , (b) iFalco et al.1 (|19991 ) 


Table 2. Observing Log 


Name 

(i) 

NED Name 
(2) 

RA 

(3) 

DEC 

(4) 

Source 

(5) 

Date 

(6) 

AM 

(7) 

texp 

(8) 

Ag 

(9) 

nyul00917 

MRK 0366 

32.8898 

13.9171 

SparsePak 

Jan 31, 2011 

1.36 

9000 

1.26 

nyu22318 

NGC 1149 

44.3495 

-0.3094 

SparsePak 

Feb 2, 2011 

1.33 

9000 

1.05 

nyu535845 

UGC 05026 

141.6659 

45.8472 

PMAS 

Jan 20, 2010 

1.08 

5400 

0.40 

nyu541044 

NGC 3921 

177.7786 

55.0788 

PMAS 

Jan 18, 2010 

1.05 

5100 

0.20 

nyu598180 

CGCG 149-037 

127.3750 

31.6755 

SparsePak 

Feb 2, 2011 

1.05 

9600 

0.68 

nyu619567 

IC 0669 

166.8190 

6.3025 

SparsePak 

Feb 1, 2011 

1.15 

9000 

0.58 

nyu674881 

NGC 6314 

258.1613 

23.2701 

CALIFA 

Apr 5, 2011 

1.09 

6300 

0.93 

nyu719486 

MRK 0385 

120.8669 

25.1027 

SparsePak 

Jan 31, 2011 

1.03 

9000 

0.44 

nyu742039 

UGC 06227 

167.8640 

47.0355 

SparsePak 

Feb 1, 2011 

1.51 

9000 

0.20 

nyu835691 

UGC 07560 

186.7364 

48.2770 

PMAS 

Jan 20, 2010 

1.02 

4600 

0.19 

nyu916578 

CGCG 216-016 

189.5908 

42.2055 

SparsePak 

Jan 31, 2011 

1.11 

9000 

0.25 

nyu916757 

NGC 4985 

197.0504 

41.6763 

SparsePak 

Jan 31, 2011 

1.02 

9900 

0.29 


Note : For each galaxy in our sample, we list (1) the ID provided in the NYU-VAGC from the SDSS DR4, (2) the most common 
associated name in the NASA Extragalactic Database, (3) the right ascension in degrees, (4) the declination in degrees, (5) the IFU on 
which that data were acquired, (6) the date of observation, (7) the average airmass of all exposures, (8) the total exposure time in 

seconds, and (9) the SDSS g-band foreground extinction in magnitudes. 


ures were taken of each of the three targets, but nyu541044 
was shorted by 5 minutes and nyu835691 by 13 minutes 
to reduce background during twilight (see Table ED. The 
second observing run was performed between January 31 st 
and February 2 nd , 2011 at the 3.5-meter WIYN telescope 
at the Kitt Peak National Observatory in Arizona us- 
ing the WIYNBench spectrograph and the SparsePak IFU 
(Bers hadv et al.l F2004I) . We used the 600 lines/mm grating 
with a blaze angle of 10.1° centred on 5500A, providing a 
spectral resolution FWHM of 5A, and wavelength range of 
3577A to 6435A. The spectra were binned 3x3, giving an in¬ 
strumental resolution of 1.4A/pixel. The 4.7" fibres provide 
a spatial resolution of 2.82 kpc at z=0.03. Five thirty-minute 
exposures were taken of each of the eight targets with ten 
and fifteen minutes additional exposure to compensate for 


reduced transparency for nyu598180 and nyu916757, re¬ 
spectively. The third run was performed as part of the Calar 
Alto Legacy Integral Field Area (CALIFA) survey on April 
11, 2011 using the Calar Alto telescope and PMAS with the 
PPAK IFU. The survey uses a three-pointing dithering with 
three exposures of 700s each per dithering with the same 
V1200 grating setup used on the previous PMAS observa¬ 
tions. In Table [2] we list relevant details from our observing 
log. 

2.3 Data Reduction 

We use the R3D reduction software (|Sanchedl2006l f to pro¬ 
duce a data cube of wavelength- and flux-calibrated spectra 
for each galaxy, which we then use to construct annular spec- 
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tra for determining the SFHs as a function of galactic radius. 
The reduction of IFU data starts with the standard proced¬ 
ures for single-slit spectra, namely subtracting the dark cur¬ 
rent and t he bias . W e pe r form cosmic ray rejection using the 
method of Ivan Dokkuml ll200lll . Because the IFU generates 
many spectra (82 for SparsePak and 382 for PMAS), the 
remaining procedures diverge from the normal methods of 
single-slit reduction. Readers interested in the technical de- 
tailsof IFU data reduction are encouraged to read ISanchej 
J2006h . Here, we provide an outline of our data reduction 
pipeline. 


2.3.1 Reduction of Fibre Data 

Locating Spectra on the Detector Each exposure from 
the IFU produces a data frame consisting of the dispersion 
axis along the horizontal, and N fibre spectra along the ver¬ 
tical (82 for SparsePak, 382 for PMAS). The location of each 
spectrum along the vertical axis is determined by finding the 
peak intensity in a sliding window of fixed size oriented ver¬ 
tically and positioned at the centre of the dispersion axis. 
Once the peaks in the centre column are identified, we next 
use them as a starting point for finding all of the other spec¬ 
tral pixels for each of the fibre spectra by recursively ap¬ 
plying the previous procedure to all of the pixels along the 
dispersion axis of the spectrograph. The fibre spectra are 
then extracted using the locations of the peak intensities 
and compiled into a row-stacked spectrum (RSS) such that 
the x-axis corresponds to the dispersion axis of the spectro¬ 
graph, and the rows along the y-axis are the fibre spectra 
listed in IFU order. 

Removing Distortions and Wavelength Calibration 

The parabolic deflection caused by the inhomogeneous dis¬ 
persion of light from the pseudoslit of the IFU fibre bundle 
as well as fibre-to-fiber distortions which are strongest on 
the edges of the detector are removed using a two-step pro¬ 
cess. First, an emission line is chosen from the arc lamp 
exposure which has strong brightness across all the fibres 
(i.e., doesn’t suffer from vignetting). The location of this 
emission line is then traced in each fibre spectrum of the 
calibration exposure, and the spectra are linearly shifted to 
a common dispersion axis such that the selected emission 
line is located at the same wavelength in each spectrum. A 
second-order correction using a fourth-order polynomial fit 
is then used to recentre all of the other emission features. 
The amount by which each spectral pixel is shifted along 
the dispersion axis is recorded in the dispersion solution. 
The observed fibre spectra are then shifted and re-entered 
in accord with the distortion solution found using the calib¬ 
ration exposure. A standard wavelength calibration is then 
applied. 

Flat-fielding, Sky Subtraction, and Flux Calibration 

A pseudospectrum is constructed for each observed fibre 
spectrum consisting of the median-combined flux at each 
wavelength from all of the fibres in the continuum exposure. 
The fibre spectra are then fibre-flattened by dividing each 
by this pseudospectrum to remove fibre-to-fiber transmis¬ 
sion, unevenness in the projection of the fibres through the 
pseudoslit, and the inhomogeneous dispersion of light from 


the pseudoslit (see ISanchej (2006) for a thorough discus¬ 
sion of these effects). The effects due to vignetting are then 
corrected by masking any pixels in the fibre-flatted spec¬ 
tra having less than seventy percent of the intensity at the 
same location in the pseudospectrum. Next, a new pseudo¬ 
spectrum is constructed by linearly interpolating the spectra 
from the sky fibres (36 for PMAS, 7 for SparsePak) onto a 
grid the same size as the IFU (382 for PMAS, and 82 for 
SparsePak), and sky-subtraction is performed on each ob¬ 
served fibre spectrum (331 for PMAS, 75 for SparsePak) us¬ 
ing these pseudospectra. Finally, the observed fibre spectra 
are flux-calibrated using spectrophotometric standard stars 
acquired during each object’s observing run. 


Constructing Data Cubes The wavelength- and flux- 
calibrated spectra are converted into a data cube by first 
populating the pixel space of the first two dimensions of 
the data cube with discretely-sampled points via a stand¬ 
ard astrometric transformation with the resulting spaxels 
(SPAtial piXELS) mapping the fibre spectra onto the sky. 
Second, a smooth, spatially-resolved spectroscopic light pro¬ 
file of each spaxel plane is constructed for each wavelength 
by blending the light of the discrete points via a two- 
dimensional Gaussian with a constant circular aperture 
(see column three of Fig. H>. Finally, galactic extinction 
along the line of sight to each target is c orrected using 
the Chariot and Fall llCharlot &; Falll l2000h law combin¬ 
ing the g-ban d extinction values given in the SDSS with 
Table 22 from Stou ghton et ah ( 2002 1 and Equation B1 from 


ISchlegel et al 


1 19981 . 


Removing Velocity Shifts In this work, we examine 
large extents of each galaxy’s surface where net rotation can 
contribute substantially to the Doppler shift. Disentangling 
the component of redshift due to rotation and recessional 
velocity is a complex task that we do not attempt here. In¬ 
stead, we fit a combination of one early-typ e and one late- 
type s pectral energy distribution (SED) from iVazdekis et ali 
(|2010f l to the spectrum at each spaxel in the data cube using 
the FIT3D code llSanchez et al.ll201lll to generate a velocity 
map containing spatially-resolved redshift information over 
the surface of the galaxy. We input into FIT3D a starting 
value of the redshift based on a spectroscopic redshift from 
either SDSS or NED (z pu bUshed from Table [T]), and the pos¬ 
sible range of allowed redshifts given by assumption of a 
maximum rotation speed of 300 km s -1 such that our meas¬ 
ured redshifts are in the range z pu bu s h e d i 300 km s 1 /c. We 
then use the resulting velocity map to remove the Doppler 
shift from the spectra. We find excellent agreement between 
our measured redshifts and those published in the SDSS or 
NED with a maximum absolute difference of ~ 30 km s _ . 
The spectra in each data cube are then shifted to a common 
rest frame such that they all begin at the same wavelength 
to ensure that co-adding spectra works correctly later. 


2.3.2 Constructing Azimuthally-binned Spectra 

In this work, we examine the stellar populations over large 
extents of each galaxy in our sample. There are several meth¬ 
ods available for combining IFU spectra to extract spatially- 
resolved stellar populations. For examp le, the CALIFA col- 
laboration use a Voronoi tessellation dCappellari fe Copinl 
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Figure 2. Images of our sample of galaxies. First column: the SDSS thumbnail images scaled to the field of view of the respective 
spectrograph from which the data were obtained with the morphological classification from M14 shown bottom centre. The white bar 
indicates 10” on the sky. Second column: the image from the first column with the binned spectral annuli overlaid and the spectrograph 
labeled. The red circle represents the half-light radius, R 50 . Third column: intensity in each spaxel of the data cube at rest frame 4050A. 
The white circle marks the largest annulus seen in the second column. 


l2003h meth od to construct local reg ions with high signal- 
to-noise ()Cid Fernandes et al.l 1 201.'ll ). We utilize a light- 
weighted azimuthal binning with constant annular width to 
provide a uniform view of the radial behaviour and con¬ 
currently counteract the decreasing signal-to-noise at large 
radii due to lower surface brightness. This method has been 
recently used to extract detailed Lyman alpha properties us¬ 
ing SSP models (jPapaderos et al1l2013l l and has been shown 
to have a nominal effect on the derived stellar population 
properties (iMast et al.l 1 2 Olfi ) . We use a two-stage process 


for measuring the S/N/pixel of our spectra. First, we use 
the standard method based on counting statistics to ensure 
only spectra with S/N above the photon-limited regime are 
considered. Importantly, we measure the S/N/pixel near the 
Balmer and 4000 break indices use d throughout this work . 
Second, we use the technique from iHusemann et al.l d2013h 
to account for the cross-correlated noise introduced by co¬ 
adding spectra. In this step, we measure the S/N/pixel in 
several continuum regions to ensure that our spectra are of 
generally high quality. 
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Error Estimate Using Counting Statistics Start¬ 
ing with the azimuthally-binned spectra, we utilize the 
signal (in counts) and the noise (in counts) at each 
wavelength to calculate the S/N/pixel. To do this, we 
first construct a row-stacked spectrum (RSS) of the es¬ 
timated Poisson noise via the standard formula noise = 
A /source + background + readout 2 where the source and 
background are the raw (i.e., not flux-calibrated) spectrum 
and sky counts, respectively, and readout is the readout 
noise from the detector. We then convert the raw spectra 
and the noise RSS into data cubes following the procedure 
outlined in Section 12.3.11 and apply the velocity shift cor¬ 
rections discussed previously. Finally, for each annulus we 
construct the azimuthal sum of the spaxels in the raw spec¬ 
tra data cube and the quadrature sum of the corresponding 
spaxels from the noise data cube to compute the S/N/pixel 
in the wavelength regions around the D n (4000), H 7 , and H<5 
indices. Because absorption features have intrinsically less 
S/N/pixel than equivalent continuum measurements, we re¬ 
tain annular spectra with a minimum S/N/pixel of 5 in each 
region around the three indices. 

Full Spectral Fitting Our analysis is centres on the equi¬ 
valent width measurements of the Balmer H 7 and H<5 in¬ 
dices. To measure these accurately, we must first determine 
the amount of nebular contamination present in our spec¬ 
tra. To measure thi s, we use the full-spectral fitting code 
FIT3D llSanchez et al.ll201lf) and a catalogue of thirty SEDs 
to find a combination which reproduces our spectra in the 
minimum chi-squared sense. We use the Single Stellar Popu¬ 
lation (SSP) S EDs generated f r om th e MILES library using 
the code from IVazdekis et al.1 ll 201 (]ll with metallicities of 
2 , 20 , and 100 percent solar each spanning a range of ten 
ages from ~ 70 Myrs to ~ 14 Gyrs. We first fit our spec¬ 
tra with all of the Balmer and forbidden oxygen emission 
lines masked. The total velocity dispersion (i.e., stellar and 
instrumental) and internal dust extinctiorfl are treated as 
free parameters and fit in the minimum chi-squared sense. 
The generated continuum fits are then subtracted from the 
spectra, and the emission lines are fit. 

We fit a single Gaussian to each emission line individu¬ 
ally with the line centre, full width half maximum, and peak 
amplitude treated as free parameters. We attempt to fit all 
of the Balmer and forbidden oxygen lines masked during the 
continuum fitting, but require that the peak of each emis¬ 
sion fit be greater than one sigma of the flux contained in a 
window of 20 around the known line centre to prevent fitting 
noise. We reject residual sky lines and any remaining cos¬ 
mic rays using a windowed boxcar filter to replace the flux 
exceeding 3 sigma of the flux in the window with the av¬ 
erage of the flux in the window. The emission fits are then 
subtracted from the original spectra, and the FIT3D con¬ 
tinuum fits are performed again with none of the emission 
lines masked. This new continuum fit is then used to per¬ 
form the emission line fitting and clipping again to produce 
a set of emission-free spectra. 

Constructing Variance Maps During the binning pro¬ 
cess spatially coherent, but physically independent spectra 

Foreground extinction has already been removed at this point. 


are co-added to increase the S/N per pixel. By combining in¬ 
dependent spectra, an artificial noise is introduced into the 
final spectrum we wish to use to make our measurements. 
To tr ace this spatially-correla ted noise, we use the meth¬ 
ods of lHusemann etahl ll2013h to construct a variance map 
to estimate the correlated noise at each wavelength using 
the SED fits to our emission-subtracted spectra calculated 
above. The variance map is simply the median-smoothed 
residuals from the full spectral fitting. Under the assump¬ 
tion that our statistical method provides the “true” meas¬ 
ure of the noise in the spectra, the variance map is con¬ 
sidered to be a direct measure of the correlated noise in¬ 
troduced during the azimuthal binning. We then use the 
variance map as the final measure of the error to recom¬ 
pute the S/N/pixel. For this final step, we use the three 
regions (4040,4060A, 4150,4200A, 4380,4430A) which meas¬ 
ure the S/N/pixel in the continuum (rather than in the 
indices) where there is less change in model-data discrep¬ 
ancies to improve our index measurements. Because we are 
measuring only continuum regions here, we require a greater 
minimum S/N/pixel than was considered when constructing 
error estimates using counting statistics where absorption 
features were used to measure the S/N/pixel. Only spectra 
with S/N/pixel > 10 are retained for further analysis. We 
show our sample with the annuli overlaid on the SDSS im¬ 
ages and the intensity map from the data cube at 4050A in 
Fig.[2]and the final set of annular spectra for each galaxy in 
Fig. ED 

2.3.3 Index Measurements 

Starting with the set of fully reduced and calibrated 
emission-subtracted annular spectra in Appendix 0 we 
measure the H<5a and H 7 A Lick indices and the narrow-band 
4000A break as a function of galactic radius. To calculate 
the strength of each index, we first define three passbands 
around each index’s line centre: the blue continuum, the 
central region, and the red continuu m. The passbands for 
the L ick H5 a and H 7 A indices from IWorthev fc Ottaviana 
(1997) and the narrow-band 4000A brea k (hereafter referred 
to as D n (4000)) from lBaloeh et all d 1999T) are shown in Fig. 
[3] To reduce random errors in the outermost annuli of our 
galaxies, we use the average Balmer index (Hc5ajH7a) = 
(H<5a + LL/a) /2 from lSanchez-Blazauez et al.[ ll2006l h 

In general, differences in D n (4000) correspond to dif¬ 
ferences in the ages of stellar populations on the order of 
several Gyrs, and differences in Balmer strength correspond 
to differing fractions of A stars. We show our index measure¬ 
ments as a function of galactic radius organized by morpho¬ 
logical classification in Fig. [4] Examining the radial trends 
of the indices, we find that strongly disturbed galaxies (visu¬ 
ally classified as SPM or pE/SPM) exhibit flat or decreasing 
Balmer absorption at increasing radii with generally increas¬ 
ing D n (4000) values with increasing galactic radius- indicat¬ 
ing that the youngest stellar populations are in the galaxies’ 
cores. For the modestly disturbed galaxies (pE) and two 
of the smooth ellipticals (E) nyu742039 and nyu619567, we 
find generally decreasing D n (4000) and increasing Balmer 
values at increasing radii consistent with their oldest stel¬ 
lar populations residing in their cores. We note that most 
of the Balmer radial indices are negative and hence not in¬ 
dicative of stellar populations containing a large fraction of 
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Figure 4. Average stellar Balmer absorption indices, (H<5 a,H7a) = (H5a + H7a) /2, and the 4000A break strength, D n (4000), as a 
function of the SDSS Petrosian r-band half-light radius (R 50 ) organized by morphological classification. Colours represent indices from 
the same galaxy in a single column. 
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Wavelength (A) 

Figure 3. Definitions of the blue continuum (upwards-hatched), 
central (horizontal-hatched), and red continuum (downwards- 
hatched) passbands for the D n (4000) (top), I[© (centre), and 
H 7 A (bottom) indices from an example galaxy. The measured 
spectrum is shown in black, and the continuum fit is shown in 
red. Note that the D n (4000) index definition does not include a 
central continuum. The y-axis shows the flux arbitrarily scaled to 
show detail, and the x-axis shows the rest frame wavelengths. 


A stars. However, there are two exceptions: nyu674881 and 
nyu916578 have positive Balmer absorption values at their 
largest radii. The final two smooth ellipticals (E) nyu598180 
and nyu535845 are distinct in that they have positive Balmer 
absorption values at nearly all radii and small D n (4000) val¬ 
ues near the range of values exhibited by the SPM galaxies. 
From Fig. [2] we see that these two galaxies possess blue cir- 
cumnuclear rings consistent with recent star formation. As 
discussed in Section [2. II our choice to include galaxies with 
interesting morphological details such as the blue rings seen 
here has paid dividends in that we find there need not be 


a direct correlation between a galaxy’s general morphology 
and its recent star formation history. 


3 MODELLING STAR FORMATION 

HISTORIES 

In this work, we seek to test if the modern merger hypo¬ 
thesis can account for the assembly of our sample of plaus¬ 
ible new ellipticals with unusually blue optical colours. Spe¬ 
cifically, we use the prediction that the blue-to-red migra¬ 
tion leaves behind a tell-tale signature via the presence of 
a strong, centrally-concentrated burst of star formation in 
the merger remnant. This signature can be detected in our 
sample of plausible new ellipticals by comparing the meas¬ 
ured Balmer and D n (4000) indices (see Section 12. 3. 3D to a 
suite of specially-constructed models. 

We use the iSEDFit tool from the impro code pack- 
agcj to construct a suite of theoretical spectral energy dis¬ 
tributions (SE Ds) covering many d ifferent SFHs (see ap¬ 
pendix A of iMoustakas et all 120131 . for details). We be¬ 
gin construction by characterizing the star formation rate 
(SFR) as a function of time for an arbitrary stellar pop¬ 
ulation using a simple, exponentially-declining SFR. such 
that SFR (t) oc e -4 / 1 " where r is the characteristic time. We 
then allow for two possible evolution paths: one in which 
the galaxy continuously forms stars throughout its lifetime 
and a second in which the galaxy experiences a burst of star 
formation superimposed atop a continuously-declining SFH 
(see Appendix 0 for details). 

The continuous models begin with a large burst of star 
formation at t = 0 and evolve for ~ 15 Gyrs with a con¬ 
tinuous, exponentially-declining SFR. For these models, we 
allow a range of e-folding times from 7 Gyrs representative 
of normal star-forming late-type galaxies to 0.5 Gyrs signi¬ 
fying a quenching of star formation. The burst models start 

https: //git hub. com / moust akas/impro 
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with the same large burst at t = 0 as in the continuous mod¬ 
els but add an additional burst of star formation generating 
between 1% and 25% of the total mass of the galaxy atop the 
exponentially-declining star formation when the galaxy is 
~ 12 Gyrs old. We use e-folding times from 25 Myrs to 1 Gyr 
representative of strongly- and weakly-truncated bursts of 
star formation, respectively. These burst models reflect the 
scenario in which our plausible new ellipticals formed from 
a recent (z = 0.12) gas-rich merger and evolved to today 
with varying degrees of quenching over a period of evolu- 
tion on the order of the dynamical merger time (~ 1 Gyr; 
iLotz et ahl l2008i l. To produce the final model spectra, we 
convolve the construct ed SFHs with a set of sin gle stellar 
population SEDs from iBruzual fe Chariot! l|2003l ) spanning 
their full range of metallicities from 0.005 Zq to 2.5 Zq with 
the Chabrier IMF. 

3.1 Index Plane 

We use the established technique of comparing the 4000A 
break strength with the Balmer Lick indices H 7 A and 
H5 a to qualitatively distinguish between galaxies that have 
a simple, continuous SFH from those that have under- 
gone a burst of star for mation within the last ~ 1 Gyr 
llKauffmann et aDl20 03a. K03). Figure [S] shows the indices 
measured from our continuous models (blue) and our burst 
models (green and red) in the Balmer-D n (4000) plane. We 
note that though there are parts of parameter space that 
can clearly only be reached by models with bursts, the burst 
models substantially overlap the regions of parameter space 
covered by the suite of continuous models; a result of the 
age-burst degeneracy. The locus of points extending along 
any continuous model track can be explained by a range of 
continuous SFHs of different e-folding times and luminosity- 
weighted ages at fixed metallicity. As shown by K03, this 
locus is not as sensitive to metallicity as the indices them¬ 
selves and acts as a simple diagnostic allowing us to qualitat¬ 
ively separate galaxies undergoing continuous star formation 
from those which have undergone a burst of star formation 
in the last 1 — 2 Gyrs in which at least 5% of their mass was 
formed. 

To simplify our analysis, we choose our fiducial burst 
and continuous models (blue, red, and green model tracks 
in Fig. E0> to have solar metallicity. We choose solar metal¬ 
licity for three reasons: 1) interpretation of the age and 
chemical composition of composite stellar populations de¬ 
termined with line-strength indices depends upon the chem¬ 
ical c omposition of the underlying ’old’ stellar popula¬ 
tions dSerra fe Tragerl 12007 1. 2) SDSS galaxies in our mass 
range have unresolved metallicities near solar dGallazzi et al.l 
12005!) . and 3) local elliptical galaxies have shallow gra dients 
i n the ir SSP-equivalent chemical compositions llTrager et al.l 
l2000al h We note, however, that within both the continuous 
and burst SFH models, metallicity gradients can plausibly 
mimic age gradients, and caution against overinterpreting 
the model plane. 

Points in the plane which lie above the continuous 
models (i.e., have a more positive Balmer absorption in¬ 
dex at a fixed D n (4000)) can be explained by a range of 
burst SFHs of different e-folding times, luminosity-weighted 
ages at fixed metallicity, and fraction of new stars formed 
during the added burst. Because of this degeneracy we do 


not try to precisely constrain model parameters such as e- 
folding times, etc. Yet, by comparing the data to our chosen 
models sampling the extrema of reasonable e-folding times 
and luminosity-weighted ages, we can draw several qualitat¬ 
ive conclusions. To do this, we divide the Balmer-D n (4000) 
parameter space into four regions by combining the known 
properties of the constructed models with the regions of the 
Balmer-D n (4000) plane where K03 found 95% of their sim¬ 
ulated galaxies had high probabilities of forming either at 
least 5% or exactly none of their stars in a burst within the 
last 2Gyrs (see figure 6 in K03 for details). 

The first region (Region 1) in the model plane con¬ 
tains the radial indices with D n (4000) < 1.35 that lie on or 
below the continuous SFH models where typical late-type 
star-forming galaxies are found. We note that points lying 
along the burst models having ages < 200 Myrs also reside in 
this region. This is an example of the well-known age-burst 
degeneracy wherein a large burst of star formation occurring 
long ago (here, the continuous model) mimics the behaviour 
of a recent, smaller burst. The second region (Region 2) 
contains the range 1.35 < D n (4000) < 1.6 and lies on or be¬ 
low the continuous SFH model with a short e-folding time 
(r = 2 Gyrs, solid blue lines). In this region, we find galaxies 
having luminosity-weighted ages between ~ 5 and 9 Gyrs, 
but note that adding a few percent of new stars by mass 
atop an older population can mimic these ages. 

The third region (Region 3) contains the radial in¬ 
dices with D n (4000) >1.6 and lying in the region where the 
continuous and burst models overlap. At fixed D n (4000), 
measurements in this region have Balmer absorption indices 
which span from the continuous SFHs consistent with old 
stellar populations in quiescent ETGs to the 5% burst mod¬ 
els with luminosity-weighted ages > 1 Gyrs. Finally, the 
fourth region (Region 4) spans D n (4000) < 1.6 and is 
bracketed above the continuous SFH models. This region 
is described by K03 as containing galaxies with a high con¬ 
fidence of having undergone a burst of star formation in the 
last ~ 2 Gyrs. The individual Balmer-D n (4000) planes with 
the radial indices for each galaxy in our sample are shown 
in Fig. IB 1 1 


3.2 Comparing Models and Data 

In Fig. H3 we overlay all of our radial index measurements 
from all galaxies in our sample atop our full suite of mod¬ 
els and find two notable conclusions. First, our measured 
indices span nearly the entire model space indicating that 
our sample of galaxies possesses a wide range of SFHs (al¬ 
though radial information is not recoverable here). Second, 
our chosen model space completely covers our measured in- 
dice4 This indicates that our choice of model parameters 
is congruent with the local population of plausible new el¬ 
lipticals. We note that many of the points that lie in the 
fourth region of the plane (i.e., bracketed above the con¬ 
tinuous models (blue)) reside below the burst models with 
burst fractions of 5% (red) and are more consistent with 
burst models with smaller burst fractions (e.g., ~ 1%). We 


The one outlier near D n (4000) ~ 1.7 is shown in Section 14.31 to 
likely be a measurement error. 
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Figure 5. The stellar Balmer absorption versus 4000A break strength model plane showing two continuous SFH models (blue), two 
burst models in which 25% of the total galaxy mass was formed in a single burst (green), and two burst models with a 5% burst fraction 
(red). Tick marks denote ages along the model lines. For the continuous models, the ages are 1, 3, 5, 7, 9, and 11 Gyrs. For the burst 
models, the ages (i.e., time since the burst) are 0.2, 0.5, 0.9, 1.5, and 2.5 Gyrs. The four qualitative regions from the text are labelled and 
orange lines mark their approximate boundaries. The full suite of continuous (top) and burst models (bottom) covering their respective 
full parameter spaces are shown in light grey. 



defer to future work to recover detailed star formation his¬ 
tories for individual data points. 

In Fig. [7] we show the Balmer and D n (4000) radial in¬ 
dices for each of our galaxies organized by the region in 
which most of the indices lie (see Appendix[B]for individual 
Balmer/D n (4000) planes). Galaxies with radial indices in 
two different regions of the plane (e.g., nyu916578) are as¬ 
signed to the region where their central annuli lie. Compar¬ 
ing with Fig. [4] we see that the first and fourth regions of 
the Balmer-D n (4000) plane contain the most morphologic¬ 
ally disturbed galaxies (T PP m ^ 10) with flat or slightly 
decreasing Balmer index gradients and generally increas¬ 
ing D n (4000) indices at large radii. With burst-like SFHs 
and very young stellar ages in their cores, these galaxies are 
likely candidates of recent gas-rich major mergers. If the dy¬ 
namical disturbance is due to merging activity as predicted 
by the modern merger hypothesis, the dichotomy of regions 
in the plane where strongly disturbed galaxies lie indicates 
that not all mergers result in star bursts- consisten t with pre¬ 
dictions from simulations (e.g., ICox et aD 12008 ). Galaxies 
in region 2 have smooth elliptical morphologies (T ppm ^ 5) 
and SFHs consistent with weak, but ongoing, star forma¬ 


tion. Galaxies in region 3 have steep positive Balmer index 
gradients and steep negative D n (4000) gradients indicative 
of very old stellar ages in their central regions, implying they 
are the least likely to be remnants of recent gas-rich major 
mergers. Surprisingly, these galaxies exhibit a variety of el¬ 
liptical morphologies (4 ^ T ppm ^ 7) including some with 
clear signs of interaction (e.g., nyu916578), but possess stel¬ 
lar ages much older than those found in the smooth ellipt¬ 
icals in region 2. 


4 DISSECTING ASSEMBLY HISTORIES 

A key prediction of the modern merger hypothesis is 
the presence of a strong, centrally-concentrated burst of 
star formation in the merger remnant. We anticipate this 
inside-out growth produces a distinct radial SFH which 
we can use to qualitatively differentiate plausible remnants 
of major mergers from galaxies undergoing other mass as¬ 
sembly mechanisms. Combining our spectroscopic radial in¬ 
dex measurements with our qualitatively defined regions of 
the Balmer-D n (4000) plane (see Section 13.11) . we ask the 
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Figure 7. The average stellar Balmer absorption indices, (TT/),\, H 74 ) = (H<5a + II 74 ) /2, and the 4000A break strength, D n (4000), 
at each annulus shown as a fraction of the SDSS Petrosian r-band half-light radius (R 50 ) for each galaxy organized by location in 
the Balmer-D n 4000 plane (see Appendix [Bj. Colours represent indices from the same galaxy in a single column. The range of T ppm 
values from Table [l] are shown for the galaxies contained in each region of the plane where larger values indicate stronger morphological 
disturbances. 


question whether each galaxy in our sample is consistent 
with the predictions of the modern merger hypothesis; i.e., 
having preferentially younger stars in its core formed in a 
recent burst. 

Additionally, we anticipate central star formation to al¬ 
ter the galaxy’s optical colours such that it either preferen¬ 
tially contains more blue light in its core or its core colours 
are reddened consistent with the presence of merger-induced 
dust. We therefore examine the distribution of global colours 
using the SDSS (u-r) and (r-z) model magnitudes (Fig. [8j 
left panel) and core colours using the SDSS fibre magnitudes 
(right panel); both are k-corrected to redshift zero in the 
AB magnitude system. For comparison, we show the control 
samples from M14 of red, spectroscopically quiescent ETGs 


and blue, star-forming LTGs from the SDSS with the same 
mass and redshift range as our sample. For the core colours, 
we show only those galaxies with z ^ 0.03 to match the red¬ 
shift range of our sample and limit the spatial bias of the 
SDSS fibre aperture]. In both panels , we show the empirical 
boundary from iHolden et alj ll2012l l which robustly separ¬ 
ates star-forming red (dusty) galaxies from non star-forming 
(old) galaxies. M14 show that 92% of all spectroscopically 
quiescent galaxies from their complete sample lie above this 
line and 97% of BPT star-formers lie below it. Additionally, 


For the redshift range of our targets, the 3” SDSS fibre aperture 
corresponds to between 0.86 and 1.8 kpc. 
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we overlay the sample of spectroscopically-confirmed star¬ 
forming ellipticals from the M14 sample as blue asterisks. 

In Fig. [31 we provide quantitative colour gradients in 
the SDSS ( u-r ) and (r-z) colour-colour plane for our galaxies 
and the subset of the comparison samples shown in Fig. [8] 
with z ^ 0.03. The colour gradients are defined by AC = 
Couter — C c ore where C is one of (u-r) or (r-z). The outer 
colour Couter is derived by converting the respective SDSS 
model magnitudes into flux, subtracting the flux from the 
SDSS fibre magnitudes, and converting to a colour. We note 
that the colour gradients computed here are not the same 
as the tracks shown in the right-hand panel of Fig. [H] which 
show the relative position of each galaxy in the colour-colour 
plane according to its core and global colours. 


4.1 Are starbursts always concurrent with major 
mergers? 


Examining Figs. [I] and [7] we find that only one (nyu541044) 
of the four galaxies visually classified as SPM or pE/SPM 
(i.e., T ppm ^ 10) has radial indices in region 4. In partic¬ 
ular, Fig. [TU] shows its central indices are consistent with 
our burst models in which ^ 15% of its stars formed in 
a recent burst and have luminosity-weighted ages ^ 1 Gyr 
(measured relative to the start of the burst). Surprisingly, 
this galaxy’s global and core optical colours both lie in the 
non star-forming region of colour-colour space as determ¬ 
ined by the Holden boundary (Fig. [8}. However, its colour 
gradients lie outside the locus of local ETGs and are more 
steeply negative (i.e., redder cores) than ~ 80% of the local 
star-forming LTGs (Fig. [9]). This combination of clear burst 
signature and quiescent core colours is consistent with post- 
starburst galaxies which are characterized as having recently 
undergone a large burst of star formation quickly followed 
by strong quenching. 

Indeed, previous spectroscopic study of this galaxy by 
ISchweizerl (j 19961 ) found this to be the case; although we 
find a slightly higher fraction of new stars formed in the 
burst (^ 15% compared to their ~ 10%) likely attribut¬ 
able to the more modern stellar template libraries used in 
constructing our burst models. Galaxies with post-starburst 
signatures are also referred to as E+A galaxies due to their 
spectra appearing as a superposition of an elliptical galaxy 
(i.e., very little nebular emission) and that of A-type stars 
lIDressler fc Gunnll 19831) . From Fig. IC1I we see the presence 
of both an elliptical-like spectrum as well as the Balmer ab¬ 
sorption features in our observed spectra for this galaxy. Ad¬ 
ditionally, the lack of dust signatures in the optical colours 
of this galaxy (Fig. [8]) is consistent with local E+A galaxies 
being quenched post-starbursts rather than starbursts with 
dust-obscured nebular emission ( lGotoll2004l ). 

The characteristics of E+A galaxies have been shown 
to be consistent with the transformation from a gas- 
rich, star-forming disk galaxy to a gas-poor , qui escent , 
pressure-supported spheroidal galaxy liNorton et al.ll200ll ). 


and are hypothesized t o come from gas-rich major mer¬ 
gers (I W ild et alj 120091 : I Cox et aT] 120081 : I Goto et aD 120081 : 
IPracv et al.l 2012|). The mol ecular gas density profile of 
this galaxy ( Hibbard fe Yunl 1 19991) is consistent with the 
centrally-concentrated gas depos ition predicted in simula - 
tions of major mergers (e.g., iBarnes fe Hernauistl fl99ll) . 
The qualitative luminosity-weighted ages we find in its cent¬ 


ral annuli are consistent with the lifetime of the short 
(~ 0.1 — 0.3 G rys) E+A signatures seen in simulations 
IlSnvder et al.11201 ll ). Without further spectroscopic meas¬ 
urements, we are unable to draw any conclusions regarding 
the role of an AGN as the truncation mechanism as pre¬ 
dicted by the modern merger hypothesis. Clearly, though, 
the preponderance of evidence indicates that this galaxy is 
a strong candidate of being a remnant of a recent, gas-rich 
major merger. 

The remaining three galaxies with T ppm ^ 10 have 
star-forming global colours with bluer core colours (Fig. 
n> and steeper positive (i.e., bluer core) colour gradients 
than ~ 96% of the local population of star-forming LTGs 
(Fig. EJ). These strongly disturbed galaxies have generally 
larger Balmer absorption values and smaller 4000A break 
strengths than all of the other morphology types (Fig. 

and specifically have their largest Balmer and smallest 
D n (4000) values in their centres- consistent with young, 
centrally-concentrated stellar populations. These character¬ 
istics match the core star formation coincident with strong 
morphological disturbances predicted by the modern mer¬ 
ger hypothesis and are consistent with these galaxies form¬ 
ing via gas-rich major mergers. Unlike nyu541044, however, 
these galaxies lack any indication of post-starburst signa¬ 
tures. Rather, they lie in region 1 (Fig. [7]) of the Balmer- 
D n (4000) plane and exhibit a dichotomy of plausible SFHs 
from long-lived continuous star formation to a small, recent 
burst of star formation. As noted in Section T3. II the overlap 
between the continuous and burst models in this region is 
a consequence of the well-known age-burst degeneracy. The 
intricacies of this effect are outside the scope of this work, 
yet we can constrain the plausible formation scenarios by 
carefully examining region 1 over our suite of models. 

In Fig. HOI we show a zoom-in of the Balmer-D n (4000) 
plane covering regions 1 and 4 with varying configurations 
of the parameter space from our models. Of the three galax¬ 
ies which lie in region 1, only nyu835691 exhibits central 
indices which are inconsistent with the continuous models. 
For the two galaxies with central annuli consistent with the 
fiducial continuous models, we note that these models ex¬ 
hibit luminosity-weighted ages in region 1 ranging between 
1 — 11 Gyrs depending on e-folding time. However, such 
large ages are in conflict with the predicted lifetimes of the 
tidal features seen in these galaxies (see below). Although 
our observations do not conclusively exclude the continuous 
models, we note that they are only plausible for these two 
galaxies and do not consider them further. Examining the 
burst models in columns 2 through 5, we see that all three 
galaxies are generally inconsistent with the extreme models 
having burst fractions of 1 and 25%, respectively. For the 
models with moderate burst fraction, we see that the ra¬ 
dial indices of all three galaxies lie along the model tracks 
with Fburst = 5% (consistent with the findings of K03 for 
this region), but only nyu719486 has central indices along 
the model tracks with the 15% burst fraction. 

Focusing on the burst models with Fburst = 5%, we find 
very young luminosity-weighted ages (< 200 Myrs, meas¬ 
ured relative to the time of the burst) over all metallicities 
and e-folding times. Given that the burst models were con¬ 
structed with the burst (and thus the merger) occurring at 
z = 0.12 and our galaxies have redshifts of z ~ 0.03, this 
leaves an intermission of ~ 1 Gyr between coalescence of 
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Figure 9. The SDSS model ( u-r ) and (r-z) colour gradients for 
massive (logM, ^ 10), nearby (z ^ 0.03) galaxies. We show the 
spectroscopically quiescent ETG (red) and blue star-forming LTG 
(blue) control groups from Fig. [8] Our sample of massive blue el¬ 
lipticals are shown using their respective morphological symbols 
and labels as in Fig. [T] Each ’blue core’ vector points in the dir¬ 
ection of its respective colour gradient having a bluer core colour. 


the merger progenitors and our observations of these galax¬ 
ies (assuming our chosen model is correct). This period is 
consistent with the lifetimes of the strong tidal signatures 
and asymmetries seen in these gala xies as w ell as in simula¬ 
tions of gas-rich major mergers jLotz et al.ll2008l ). 

The lack of clear signatures of a large burst of star form¬ 
ation in three of the four SPM or pE/SPM (i.e., T ppm ^ 10) 
galaxies in our sampl e is consistent with previous studies of 
blue spheroids fe.g.. iToieiro et aL |2013l ) and actively star¬ 
forming local ellipticals (Fukueit a et al.l 120041 ) which were 
shown to have SFHs similar to the long-lived star formation 
seen in local star-forming LTGs. This indicates that not all 
plausible major merger remnants undergo an intense burst of 
star formation. Indeed, this dichotomy is also seen in simula¬ 
tions of gas-rich major mergers which show that the strength 
and duration of merger-induced starburst activity is likely 
dependent upon the gas mass fra ction and mass ratios of the 
progenitor disks (I Cox et al . 2008) as well as orbital paramet¬ 
ers at fixed masses ( Snyder et al.l 120111 ). We therefore find 
that they are likely candidates of being remnants of recent, 
gas-rich major mergers. 


4.2 Dusty ellipticals temporarily removed from 
the red sequence? 

As shown in Fig. [7] galaxies with radial indices lying princip¬ 
ally in region 2 of the Balmer-D n (4000) plane exhibit relat¬ 
ively flat gradients in D n (4000) (AD n (4000) ~ 0.1) but non- 
uniform gradients in Balmer absorption that peak at radii 
> lR e . The flat gradients in D n (4000) with variable Balmer 
absorption are indicative of a changing fraction of young 
(< 1 Gyrs) stars atop a uniform, older popu l ation . Such a 
“frosting” of low-level star format ion (ITrager et al.l l2000bl : 
lYi et al. 120051 : iKavirai et al.ll2009l : ISalim fc Rich 20101) may 
add only 1—2% to a galaxy’s mass, but has enough hot O and 
B stars to produce a brief blue colour migrating these ellipt¬ 
icals temporar ily blue-ward off of the red sequence into th e 
green valley (ICortese fe Hughesl 120091 : iThilker et al] l20ld) . 
Overall, these galaxies have smooth morphologies (E), but 
nyu674881 is labelled pE owing to its large, circumgalactic 
dust ring (Fig. 0. 


Examining their colours in Fig. [8] we find the two ellipt¬ 
icals (T ppm = 5) have global colours which lie in the region 
of the colour-colour plane occupied by star-forming LTGs 
as well as the spectroscopically-confirmed star-forming el¬ 
lipticals from M14. Their core colours track upward along 
the dust-reddening vector away from the star-forming el¬ 
lipticals toward redder colours; suggesting the presence of 
central dust obscuration. Their colour gradients (Fig. [9j lie 
outside the locus of quiescent ETGs with (u-r) gradients 
a factor of ~ 2 steeper than those seen in the quiescent 
ETGs. Such steepened colour gradients (i.e., redder core 
colours) have been previously identified in nearby, lu min¬ 
ous elliptical galaxies with optical dust features (iKirn fe Iml 
120131) and are likel y seede d from accretion of a gas-rich satel¬ 
lite (iMartini et al.l 120131 ). Additionally, these galaxies pos¬ 
sess circumnuclear blue rings (Fig. [2j associated wit h rej u- 
venation of quiescent stellar systems IlShapiro et al.l 120101) . 
consistent with a frosting interpretation. 

Both the global and core colours of the dust-ring pE 
galaxy (nyu674881) reside in the gap between star-forming 
and quiescent galaxies. Its relatively flat colour gradient is 
consistent with local quiescent ETGs; the dust ring appears 
to have little effect on the its colour gradient in contrast 
to the other two galaxies in this region. Yet, previous ana¬ 
lysis of elliptical galaxies with prominent optical dust lanes 
(rings) indicates that these features represent a transition 
phase between starburst and quiescence and implicate re¬ 
cent gas-rich minor (i.e., Mi,« : M 2 .« > 4) merging as their 
likely cause llShabala et al.[|2012l ). The capture of a gas-rich 
satellite via minor merging has been shown to induce low- 
level star formatio n cons istent with a frosting interpreta¬ 
tion l|Kavirai et al.ll2009l ). We note that if minor merging 
is responsible for the enhanced star formation, the galaxy’s 
generally smooth morphology indicates that the star form¬ 
ation lifetime appears to be greater than the time for the 
interaction features to disappear. Overall, all three galaxies 
from region 2 of the Balmer-D n (4000) plane lack the clearly 
young spectral signatures and blue core colours seen in the 
more disturbed SPM and pE/SPM galaxies and are hence 
not as likely to be remnants of recent major mergers. 


4.3 Many Paths into the Green Valley? 

The five galaxies that principally lie in the third region 
of the Balmer-Dn (4000) plane exhibit a tight clustering of 
both radial indices with a generally negative (i.e., smal¬ 
ler values at larger radii) slope in their 4000A break 
strengths and a positive slope in their Balmer absorp¬ 
tion values (Fig. 0. The central D n (4000) index values in 
these galaxies are the largest seen at any radius of any 
galaxy in our sample suggesting that these galaxies are 
the least likely to be remnants of recent gas-rich major 
mergers due to old (> 5 Gyrs) populations in their cores. 
The overall decreasing 4000A break strength and increas¬ 
ing (but still negative) Balmer absorption at large radii 
are consistent with population gradients seen i n the genera l 
population of nearby E TGs IlGonzalez Delgado et al. [2014; 

ISanchez-Blazquez et al.l 120061) . The anomalous large, posit¬ 
ive Balmer absorption value seen at the sixth annulus of 
nyu619567 is likely a measurement error stemming from an 
over-fitting of the H5 a absorption feature (see Fig. ICTl) . 

These galaxies have quiescent global and core colours 
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Figure 8. [left] Global ( u-r ) and ( r-z ) colours and [right] core colours from the SDSS k-corrected to redshift zero in the AB magnitude 
system over-plotted on the quiescent (red) and star-forming (blue) control samples described in Section [4] Core colours are shown only 
for galaxies with z ^ 0.03 to match the redshift of our sample and minimize aperture effects. The dotted lines in the figure on the right 
connect the core colours with the correspond ing global colours for each galaxy. Symbols represent visual classifications from M14. The 
empirical boundary from Hol den et al.l d20121) (green) robustly separates the star-forming and non star-forming galaxies with star-forming 
galaxies below the line. The complete sample of star-forming elliptical galaxies from M14 is shown in blue asterisks. Letter labels are 
from Table [I] 


with small, negative gradients (i.e., slightly redder core 
colours, Fig. 0 consistent with the locus of local qui¬ 
escent ETGs. Their stellar massed are in excess of 3 x 
10 10 Mq- the characteristic mass above which local galax¬ 
ies tend to have little star forma tion and increasingly be - 
come bulge-dominated ellipticals (IKauffmann et alJl2003bli . 
These galaxies have a median scatter in their ( g-r ) colours 
blue-ward from the empirical red/blue dividing line in Fig. 
[Q of ~ 0.03 mags. This is within the typical colour error 
(0.04 mags) found by M14 and may account for their appar¬ 
ent green-valley colours. It is also possible that these galax¬ 
ies, like those found in region 2, may be red sequence galaxies 
which have migrated into the green valley. Notably, however, 
they lack the Balmer absorption seen in the region 2 galaxies 
and have luminosity-weighted ages at least 1 — 2 Gyrs older. 

The galaxies in region 3 present a rather homogeneous 
group reminiscent of local relaxed ellipticals with smooth 
cores, quiescent colours, steep index gradients, and no dust 
features. Further, all of the elliptical (i.e., not SPM) galaxies 
in our sample classified as pE owing to outer asymmetries 
(e.g., nyu916578) lie in this region. This is surprising as we 
argued previously that the galaxies in region 2 are likely 
to be remnants of minor mergers, but those galaxies lack 
the structural signatures of recent interaction (i.e., arms, 
loops, shells, etc.) seen in these pE galaxies. This raises the 
important question as to how the galaxies in region 3 could 
plausibly come about. 

It is well known that mergers between equal-mass disk 
galaxies result in remnants with tidal tails and plumes 
llToomre fc Toomr j 1 19721 : iBarnesI 1 19881 . 1 199211 , but the low 
surface brightness of these and other late-stage features (e.g., 
small loops, arms, and shells) are also seen in minor mergers 
llFeldmann et al.ll2008h making discriminating b etween these 
two types of merger remnants difficult l|Lotz et al.l 120081 : 
IStruck fe Smithl l2012h . Further, not all morphologically 


Using h = 0.7. 


distu rbed ellipticals are star-forming (jMichard fe PrugnieJ 
l2004h , and not all blue (presumably star-forming) ellipticals 
are morphologically disturbed (see nyu22318 in our sample). 
Yet both plausible formation histories predict the low-level 
star formation observed in nearby ETGs (e.g., lYi et al] 
12005h . Additionally, a scenario in which an existing ellipt¬ 
ical galaxy undergoes a (nearly) equal-mass gas-rich mer¬ 
ger would leave behind a remnant with such late-stage fea¬ 
tures. These so-called “mixed” mergers have been shown to 
be impor tant in explaining the dich otomy of local elliptical 
galaxies ijRothberg fc Josephl l2006bl ). However, their gen¬ 
eral imp ortance in g alaxy evolution is not well constrained 
llKhochfar fc Burkertll2003l l. 

The modestly disturbed morphologies and quiescent 
colours of the galaxies in region 3 may be described by 
yet another merger scenario: dissipationless (so-called “dry”) 
major mergers. In this scenario, two gas-poor elliptical 
galaxies merge, producin g a rem n ant wi th broad, low sur¬ 
face brightness features (IBell et al.l 12003 ) like those seen in 
many of the pE galaxies in our sample. Assuming the pro¬ 
genitor galaxies posses colours slightly redder than our em¬ 
pirical red/blue dividing line, such a dry merger would con¬ 
serve colour but allow the remnant to cross our dividing line; 
producing a slightly bl uer-than- usual elliptical galaxy at a 
given stellar mas s (ISkelt on et al.ll2009l 'l However, simulations 
from iKawata et all (l2006i l show that the colour and struc¬ 
ture of such features can also be produced through spheroid- 
disk gas-rich minor mergers. The same complications lis¬ 
ted above for m i xed-mergers th en a pply to this scenario, as 
well. ISanchez-Blazq ugz_et_al.l j|2009|) found that half of the 
ellipticals from van Dokkuml (120051 1 having strong tidal fea¬ 
tures exhibit a “frosting” of star formation like that seen in 
our sample of galaxies in region 2. These frosted ellipticals 
were also found to be supported by rotation, indicating the 
presence of a dynamically cold stellar component. This is in 
contrast to simulations of dry mergers which produce rem¬ 
nants with dispersion-do minated kinematics similar to those 
of local cluster galaxies dNaab et al.[i2006l ; iNaab fc Ostrikeil 
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Figure 10. Zoom-in of regions 1 and 4 of the Balmer-D n (4000). Column 1 shows the full suite of continuous models in grey and highlights 
the two fiducial models used throughout this work in blue. Columns 2 through 5 show the full suite of burst models in grey with the 
fraction of stellar mass formed during the burst increasing to the right. The two sets of fiducial burst models used throughout this work 
are shown in red (5% burst) and green (25% burst). The fiducial continuous model with an e-folding time of 7Gyrs is shown in solid 
dark grey alongside the burst models as a guide. Each row shows the radial index measurements corresponding to the galaxy listed in 
the first column. 



120091) . While our sample of galaxies in region 3 is incon¬ 
sistent with the frosting scenario, it follows that a careful 
analysis of the kinematic structure of these galaxies could 
be used to plausibly identify if dry mergers may play a role 
in their assembly. 

It is clear, then, that there are many routes a galaxy 
may take to arrive in region 3 of the Balmer-D n (4000) plane. 
Extracting the precise star formation history for these galax¬ 
ies requires analysis outside the scope of this work. However, 
we note that in Fig. Q the galaxy nyu916578 stands out 
among the others in region 3 in that its three outermost 
radial indices track into the fourth region of the Balmer- 
Dn (4000) plane near the burst models with luminosity- 
weighted ages of ~ 1 Gyr and exhibit Balmer absorption 
consistent with an increased fraction of young stars. This 
is highly contrasted by its core population with ages closer 


to ~ 10 Gyrs. Such population gradients have been seen 
in “revitalized” elliptical galaxies consistent with ongoing 
star formation in their outer annuli ass ociated with tida l 
stripping during a gas-rich minor merger dFang et al.ll2012l ). 
The presence of younger stellar populations at large galactic 
radii is consistent with the outside-in formation scenario 
of elliptical galaxies and offers a possible explanation for 
observed "inverted" age gradients and lower m etallicity at 
large radii in some local ellip tical galaxies (e.g., iBaes et~akl 
120071 : lLa Barbera et akll20ld ) . 


5 DISCUSSION 

Using the predictions of the modern merger hypothesis as a 
guide, we measure the radial Balmer and D n (4000) indices 
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to attempt to distinguish between galaxies originating from 
a recent, gas-rich major merger and those undergoing some 
other mass assembly mechanism. In Section 3] we present a 
detailed analysis combining the radial indices, optical col¬ 
ours, and morphological features to examine the star form¬ 
ation histories of our sample of galaxies. With those details 
in mind, we show a characteristic “decision” tree (Fig. 1111) 
for our sample which allows us to a posteriori correlate the 
characteristics of each galaxy in our sample with its location 
in the qualitatively defined regions of the Balmer-D n (4000) 
plane. 

We begin by splitting our sample into two groups us¬ 
ing the presence of strong outer loops. This cleanly separates 
strongly disturbed spheroidal post merger (SPM) signatures 
(i.e., T ppm ^ 10) from ellipticals with smooth or modestly 
disturbed morphologies. For the most disturbed galaxies, we 
further examine whether they are consistent with undergo¬ 
ing a strong central burst of star formation. As discussed in 
Section mu the key distinction between galaxies residing in 
regions 1 and 4 appears to be only the strength of the burst 
and not the existence of a burst. Further, we find that galax¬ 
ies in region 4 appear to be consistent with E+A or post- 
starburst galaxies. However, we note that our small sample 
statistics preclude extending this result further. Galaxies in 
region 1 of the index plane are consistent with a dichotomy 
of star formation histories from long-lived, continuous star 
formation to a recent, weak burst, forming ^ 5% of their 
mass in a recent burst. 

For the second half of the tree, we find that galaxies 
with smooth elliptical morphologies but internal structures 
such as circumnuclear blue rings or dust lanes preferentially 
reside in region 2 of the index plane. Additionally, they pos¬ 
sess flat gradients in their 4000A break strengths, but have 
radial Balmer absorption indices consistent with a changing 
fraction of young stars atop a uniform, older population. 
As outlined in Section mi previous examination of galax¬ 
ies with similar internal structures attribute these features 
to revitalized stellar populations brought on by minor mer¬ 
ger events- in agreement with our frosting interpretation. 
This is in accord with the general idea that minor mer¬ 
gers are thought to play a significant role in evolving spher¬ 
oidal galaxies for two reasons: (i) the small number of ma¬ 
jor mergers in the local uni verse cannot prov i de the abund- 
ance of green-valley E TGs (iDarg et al.l l 2 Qlol : ISalim fe Richl 
120101: iLotz et al.ll201lh . and (ii) the observed size evolution 
of elliptical galaxies since z ~ 1 i s best ex plai ned b y growth 
dominated by minor mergers llBedorf fc Porteeies Zward 
2013^ Bezanson_etalJ 2009^ Hilzet ahl 20131 : iLee et al.ll2013l : 


Naab et, alJl 200 fll : lOser et, al.ll 2012 l i. 


Moving to the final branch of the decision tree, we an¬ 
ticipate finding smooth (i.e., small T ppm values) galaxies 
with star formation histories least like those predicted in 
the modern merger hypothesis. Indeed, all but one of the 
galaxies in region 3 have star formation histories much like 
local quiescent ETGs. However, many of them exhibit ex¬ 
ternal structures such as arms or loops, and all of the ellipt¬ 
ical (i.e., not SPM) galaxies labelled as pE owing to these 
outer asymmetries (e.g., nyu916578) lie in this region. The 
one galaxy which runs counter to this general interpret¬ 
ation is nyu916578 which exhibits a strong, negative age 
gradient with stellar populations at its largest radii having 
luminosity-weighted ages < 1 Gyr. Such “inverted” age gradi¬ 


ents are attributable to gas-rich minor merging. Overall, the 
galaxies in region 3 have a scatter in their ( g-r ) colours of 
only ~ 0.03 mags blue-ward of the empirical red/blue divid¬ 
ing line (Fig. [lj and lack clear indications of ongoing star 
formation. It is likely that their green valley colours will 
quickly fade and migrate them onto the red sequence. Ad¬ 
ditionally, they possess stellar masses near ~ 10 10 ' 8 Mq1i~ 2 
and extinction- and k-corrected SDSS absolute r-band mag¬ 
nitudes between —21 ^ M r ^ —22 (Table Q, placing them 
close to the knee of the green valley luminosity function 
found bv lGoncalves et all 1 2012|j . consistent with the downs¬ 
izing of the red sequence (i.e., the low-mass end of the red 
sequence is populated at late times). 

Previous detailed examination of green valley 
galaxies revealed a d ichotomy in their morphologies 
IlSchawinski et al.l l2014'). The presence of both early- and 
late-type galaxies in the green valley lends support to 
the idea that multiple mass-assembly mechanisms are 
likely responsible for galaxy growth over cosmic time. Our 
targeted study of spheroidal galaxies with blue sequence 
colours expands on that work, demonstrating a further 
branching of mass assembly histories even within just 
early-type galaxies in the green valley. This indicates 
that the green valley is a rich laboratory for studying 
the transitory phases of galaxy formation and evolution. 
Our analysis could only be accomplished via the usage of 
spatially-resolved spectroscopy and motivates the study of 
green valley galaxies in large IFU surveys. 


6 SUMMARY AND CONCLUSIONS 

In this work, we have combined visual morphologies with 
spatially-resolved spectra and optical colours to character¬ 
ize the radial SFHs of a unique sample of plausible new el¬ 
lipticals. Utilizing techniques for examining SFHs previously 
used to characterize the central regions of local galaxies (e.g., 
K03), we compare the stellar age indices H<5, H 7 , and the 
4000A break to synthetic SFH models at radii from the core 
out to ~ 2 — 3 half-light radii. We find that these index 
values and their radial dependence correlate with specific 
morphological features such that the most disturbed galax¬ 
ies have the smallest 4000A break strengths and the largest 
Balmer absorption values. With this new spatial data, we 
are able to detect the presence of a centrally concentrated 
starburst occurring within the last ~ 1 Gyr - allowing us 
to plausibly discriminate between the construction of new 
ellipticals by gas-rich major merging following the modern 
merger hypothesis (Hopkins ct al. 12008a) and ellipticals un¬ 
dergoing other mass assembly mechanisms. 

Our main conclusions are as follows. 

(i) We find that one-third of our sample possess the strong 
morphological disturbances and central star formation con¬ 
sistent with being plausible remnants of a recent, gas-rich 
major merger. However, only one of them exhibits the char¬ 
acteristics of having recently undergone a strong, centrally- 
concentrated burst of star formation followed by quench¬ 
ing. Its quiescent optical colours and strong burst signature 
are consistent with its well-known status as an E+A (post- 
starburst) galaxy. Assuming our small sample comports with 
the population of local merger remnants, we infer that the 
modern merger hypothesis describes the assembly of only 
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Figure 11. A characteristic “decision” tree for the classification of our sample of massive, blue spheroidal galaxies. The images come from 
the SDSS and are scaled to 40kpc/h on a side. At the top of each image, we list the name, single-letter designation, and morphological 
classification from Table |T] The T ppm and morphological characteristics are detailed in Section 12.11 Each large box represents one of the 
regions of the Balmer-D n (4000) plane as described in Section 13.II 


~ 25% of plausible major merger remnants in the local uni¬ 
verse. 

(ii) Three of the four plausible major mergers in our 
sample have star-forming optical colours with much bluer 
centres. Their Balmer absorption and 4000A break strength 
index values indicative of having centrally-concentrated 
young stellar populations. However, they lack the clear sig¬ 
natures of a recent, strong central burst of star formation. 
Overall, their star formation histories are generally consist¬ 
ent with having formed ~ 5% of their stellar mass in a re¬ 
cent, weak burst. 

(iii) Two-thirds of our sample of 12 massive (M* ^ 
10 10 Mq), nearby (z ^ 0.03), visually-selected plausible new 
ellipticals with blue-sequence optical colours and varying de¬ 
grees of morphological peculiarities are inconsistent with the 
predictions from simulations of gas-rich major mergers. We 
therefore find that their mass assembly histories are not well- 
described by the modern merger hypothesis. Half of these 
galaxies exhibit characteristics seen in remnants of recent, 


gas-rich minor mergers and the other half have radial star 
formation histories consistent with those seen in local qui¬ 
escent ellipticals. 

We are happy to thank Hans-Walter Rix, Anna Gallazzi, 
and Enrique Perez for useful insights and discussions. T.H. 
acknowledges support from the Missouri Consortium of 
NASA’s National Space Grant College and Fellowship Pro¬ 
gram. D.H.M. acknowledges support from the Research 
Corporation for Science Advancement under the Cottrell 
College Science Award grant No. 10777. Funding for the 
SDSS has been provided by the Alfred P. Sloan Founda¬ 
tion, the Participating Institutions, the National Aeronaut¬ 
ics and Space Administration, the National Science Found¬ 
ation, the U.S. Department of Energy, the Japanese Mon- 
bukagakusho, and the Max Planck Society. The SDSS Web 
site is http://www.sdss.org/ The SDSS is managed by the 
Astrophysical Research Consortium (ARC) for the Parti¬ 
cipating Institutions, which are The University of Chicago, 
Fermilab, the Institute for Advanced Study, the Japan Parti- 


© 0000 RAS, MNRAS 000, 000-000 









18 T. Haines et al. 


cipation Group, The Johns Hopkins University, Los Alamos 
National Laboratory, the Max-Planck-Institute for Astro¬ 
nomy (MPIA), the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, University of Pitts¬ 
burgh, Princeton University, the United States Naval Ob¬ 
servatory, and the University of Washington. This research 
is (partially) based on data from the MILES project. This 
publication also made use of NASA’s Astrophysics Data Sys¬ 
tem Bibliographic Services. Based on observations collec¬ 
ted at the Centro Astronomico Hispano Aleman (CAHA) at 
Calar Alto, operated jointly by the Max-Planck Institut fur 
Astronomie and the Instituto de Astrofisica de Andalucia 
(CSIC). This paper is (partially) based on data obtained 
by the CALIFA Survey, funded by the Spanish Ministry of 
Science under grant ICTS-2009-10 [and add-ons], and the 
Centro Astronomico Hispano-Aleman. 


References 

Adelman-McCarthy J. K., et al., 2006, The Astrophysical 
Journal Supplemental, 162, 38 
Bacon R., et al., 2001, Monthly Notices of the Royal As¬ 
tronomical Society, 326, 23 

Baes M., Sil’chenko O. K., Moiseev A. V., Manakova E. A., 
2007, Astronomy and Astrophysics, 467, 991 
Baldry I. K., Glazebrook K., Brinkmann .J., Ivezic v., 
Lupton R. H., Nichol R. C., Szalay A. S., 2004, The As¬ 
trophysical Journal, 600, 681 

Balogh M. L., Morris S. L., Yee H. K. C., Carlberg R. G., 
Ellingson E., 1999, The Astrophysical Journal, 527, 54 
Barnes J. E., 1988, The Astrophysical Journal, 331, 699 
Barnes J. E., 1992, The Astrophysical Journal, 393, 484 
Barnes J. E., Hernquist L. E., 1991, The Astrophysical 
Journal Letters, 370, L65 

Barnes J. E., Hernquist L., 1996, The Astrophysical 
Journal, 471, 115 

Bedorf J., Portegies Zwart S., 2013, Monthly Notices of the 
Royal Astronomical Society, 431, 767 
Bell E. F., McIntosh D. H., Katz N., Weinberg M. D., 2003, 
The Astrophysical Journal Supplement Series, 149, 289 
Bell E. F., et al., 2004a, The Astrophysical Journal Letters, 
600, Lll 

Bell E. F., et al., 2004b, The Astrophysical Journal, 608, 
752 

Bell E. F., et al., 2006, The Astrophysical Journal, 640, 241 
Bell E. F., Zheng X. Z., Papovich C., Borch A., Wolf C., 
Meisenheimer K., 2007, The Astrophysical Journal, 663, 
834 

Bell E. F., et al., 2012, The Astrophysical Journal, 753, 167 
Bershady M. A., Andersen D. R., Harker J., Ramsey L. W., 
Verheijen M. A. W., 2004, Publications of the Astronom¬ 
ical Society of the Pacific, 116, 565 
Bezanson R., van Dokkum P. G., Tal T., Marchesini D., 
Kriek M., Franx M., Coppi P., 2009, The Astrophysical 
Journal, 697, 1290 

Blanton M. R., 2006, The Astrophysical Journal, 648, 268 
Blanton M. R., et al., 2005, The Astronomical Journal, 129, 
2562 

Brammer G. B., et al., 2009, The Astrophysical Journal 
Letters, 706, L173 

Brown M. J. I., Dey A., Jannuzi B. T., Brand K., Benson 


A. J., Brodwin M., Croton D. J., Eisenhardt P. R., 2007, 
The Astrophysical Journal, 654, 858 
Bruzual G., Chariot S., 2003, Monthly Notices of the Royal 
Astronomical Society, 344, 1000 
Cappellari M., Copin Y., 2003, Monthly Notices of the 
Royal Astronomical Society, 342, 345 
Cappellari M., et al., 2011, Monthly Notices of the Royal 
Astronomical Society, 413, 813 
Chariot S., Fall S. M., 2000, The Astrophysical Journal, 
539, 718 

Cid Fernandes R., et al., 2013, Astronomy & Astrophysics, 
557, A86 

Cortese L., Hughes T. M., 2009, Monthly Notices of the 
Royal Astronomical Society, 400, 1225 
Cox T. J., Jonsson P., Primack J. R., Somerville R. S., 
2006, Monthly Notices of the Royal Astronomical Society, 
373, 1013 

Cox T. J., Jonsson P., Somerville R. S., Primack J. R., 
Dekel A., 2008, Monthly Notices of the Royal Astronom¬ 
ical Society, 384, 386 

Darg D. W., et al., 2010, Monthly Notices of the Royal 
Astronomical Society, 401, 1043 
Dressier A., Gunn J. E., 1983, The Astrophysical Journal, 
270, 7 

Eliche-Moral M. C., et al., 2010, Astronomy and Astro¬ 
physics, 519, 55 

Faber S. M., et al., 2007, The Astrophysical Journal, 665, 
265 

Falco E. E., et al., 1999, Publications of the Astronomical 
Society of the Pacific, 111, 438 
Fang J. J., Faber S. M., Salim S., Graves G. J., Rich R. M., 
2012, The Astrophysical Journal, 761, 23 
Feldmann R., Mayer L., Carollo C. M., 2008, The Astro- 
physical Journal, 684, 1062 

Fukugita M., Nakamura O., Turner E. L., Helmboldt J., 
Nichol R. C., 2004, The Astrophysical Journal Letters, 
601, L127 

Gallazzi A., Chariot S., Brinchmann J., White S. D. M., 
Tremonti C. A., 2005, Monthly Notices of the Royal As¬ 
tronomical Society, 362, 41 

Gon§alves T. S., Martin D. C., Menendez-Delmestre K., 
Wyder T. K., Koekemoer A., 2012, The Astrophysical 
Journal, 759, 67 

Gonzalez Delgado R. M., et al., 2014, Astronomy & Astro¬ 
physics, 562, A47 

Goto T., 2004, Astronomy and Astrophysics, 427, 125 
Goto T., Kawai A., Shimono A., Sugai H., Yagi M., Hat- 
tori T., 2008, Monthly Notices of the Royal Astronomical 
Society, 386, 1355 

Hibbard J. E., Yun M. S., 1999, The Astrophysical Journal 
Letters, 522, L93 

Hilz M., Naab T., Ostriker J. P., 2013, Monthly Notices of 
the Royal Astronomical Society, 429, 2924 
Holden B. P., van der Wei A., Rix H.-W., Franx M., 2012, 
The Astrophysical Journal, 749, 96 
Hopkins P. F., Hernquist L., Cox T. J., Robertson B., 
Springel V., 2006, The Astrophysical Journal Supplement 
Series, 163, 50 

Hopkins P. F., Hernquist L., Cox T. J., Keres D., 2008a, 
The Astrophysical Journal Supplement Series, 175, 356 
Hopkins P. F., Cox T. J., Keres D., Hernquist L., 2008b, 
The Astrophysical Journal Supplement Series, 175, 390 


© 0000 RAS, MNRAS 000, 000-000 



Testing the Modern Merger Hypothesis 19 


Husemann B., et al., 2013, Astronomy and Astrophysics, 
549, 87 

Kauffmann G., White S. D. M., Guiderdoni B., 1993, 
Monthly Notices of the Royal Astronomical Society, 264, 
201 

Kauffmann G., et ah, 2003a, Monthly Notices of the Royal 
Astronomical Society, 341, 33 
Kauffmann G., et al., 2003b, Monthly Notices of the Royal 
Astronomical Society, 341, 54 
Kaviraj S., Peirani S., Kliochfar S., Silk J., Kay S., 2009, 
Monthly Notices of the Royal Astronomical Society, 394, 
1713 

Kawata D., Mulchaey J. S., Gibson B. K., Sanchez- 
Blazquez P., 2006, The Astrophysical Journal, 648, 969 
Kelz A., et al., 2006, Publications of the Astronomical So¬ 
ciety of the Pacific, 118, 129 

Kliochfar S., Burkert A., 2003, The Astrophysical Journal 
Letters, 597, L117 

Kliochfar S., Silk J., 2006, The Astrophysical Journal Let¬ 
ters, 648, L21 

Kim D., Im M., 2013, The Astrophysical Journal 
Kormendy J., 2009, in The Elliptical-Spheroidal and 
Elliptical-Elliptical Galaxy Dichotomies, p. 87 
La Barbera F., De Carvalho R. R., De La Rosa I. G., Gal 
R. R., Swindle R., Lopes P. A. A., 2010, The Astronomical 
Journal, 140, 1528 

Lee J. H., Lee M. G., Hwang H. S., 2006, The Astrophysical 
Journal, 650, 148 

Lee J. H., et al., 2013, The Astrophysical Journal Letters, 
762, L4 

Lotz J. M., Jonsson P., Cox T. J., Primack J. R., 2008, 
Monthly Notices of the Royal Astronomical Society, 391, 
1137 

Lotz J. M., Jonsson P., Cox T. J., Croton D., Primack 
J. R., Somerville R. S., Stewart K., 2011, The Astrophys¬ 
ical Journal, 742, 103 

Martin D. C., et al., 2007, The Astrophysical Journal Sup¬ 
plement Series, 173, 342 

Martini P., Dicken D., Storchi-Bergmann T., 2013, The As¬ 
trophysical Journal, 766, 121 

Mast D., et al., 2014, Astronomy & Astrophysics, 561, A129 
McIntosh D. H., et al., 2014, Monthly Notices of the Royal 
Astronomical Society, 442, 533 
Michard R., Prugniel P., 2004, Astronomy and Astrophys¬ 
ics, 423, 833 

Mihos J. C., Hernquist L., 1996, The Astrophysical 
Journal, 464, 641 

Moustakas J., et al., 2013, The Astrophysical Journal, 767, 
50 

Naab T., Ostriker J. P., 2009, The Astrophysical Journal, 
690, 1452 

Naab T., Kliochfar S., Burkert A., 2006, The Astrophysical 
Journal Letters, 636, L81 

Naab T., Johansson P. H., Ostriker J. P., 2009, The Astro- 
physical Journal Letters, 699, L178 
Norton S. A., Gebhardt K., Zabludoff A. I., Zaritsky D., 
2001, The Astrophysical Journal, 557, 150 
Oser L., Naab T., Ostriker J. P., Johansson P. H., 2012, 
The Astrophysical Journal, 744, 63 
Papaderos P., et al., 2013, Astronomy and Astrophysics, 
555, LI 

Pracy M. B., Owers M. S., Couch W. J., Kuntschner H., 


Bckki K., Briggs F., Lah P., Zwaan M., 2012, Monthly 
Notices of the Royal Astronomical Society, 420, 2232 
Rosales-Ortega F. F., Kennicutt R. C., Sanchez S. F., 
Diaz A. I., Pasquali A., Johnson B. D., Hao C. N., 2010, 
Monthly Notices of the Royal Astronomical Society, 405, 
735 

Roth M. M., et al., 2005, Publications of the Astronomical 
Society of the Pacific, 117, 620 
Rothberg B., Joseph R. D., 2006a, The Astronomical 
Journal, 131, 185 

Rothberg B., Joseph R. D., 2006b, The Astronomical 
Journal, 132, 976 

Salim S., Rich R. M., 2010, The Astrophysical Journal Let¬ 
ters, 714, L290 

Sanchez S. F., 2006, Astronomische Nachrichten, 327, 850 
Sanchez-Blazquez P., Gorgas J., Cardiel N., 2006, Astro¬ 
nomy and Astrophysics, 457, 823 
Sanchez-Blazquez P., Gibson B. K., Kawata D., Cardiel 
N., Balcells M., 2009, Monthly Notices of the Royal As¬ 
tronomical Society, 400, 1264 
Sanchez S. F., et al., 2004, The Astrophysical Journal, 614, 
586 

Sanchez S. F., Rosales-Ortega F. F., Kennicutt R. C., John¬ 
son B. D., Diaz A. I., Pasquali A., Hao C. N., 2011, 
Monthly Notices of the Royal Astronomical Society, 410, 
313 

Sanchez S. F., et al., 2012, Astronomy and Astrophysics, 
538, 8 

Schawinski K., Thomas D., Sarzi M., Maraston C., Kaviraj 
S., Joo S.-J., Yi S. K., Silk J., 2007, Monthly Notices of 
the Royal Astronomical Society, 382, 1415 
Schawinski K., et al., 2009, Monthly Notices of the Royal 
Astronomical Society, 396, 818 
Schawinski K., Dowlin N., Thomas D., Urry C. M., Ed¬ 
mondson E., 2010, The Astrophysical Journal Letters, 
714, L108 

Schawinski K., et al., 2014, Monthly Notices of the Royal 
Astronomical Society, 440, 889 
Schlegel D. J., Finkbeiner D. P., Davis M., 1998, The As¬ 
trophysical Journal, 500, 525 

Schweizer F., 1996, The Astronomical Journal, 111, 109 
Serra P., Trager S. C., 2007, Monthly Notices of the Royal 
Astronomical Society, 374, 769 
Shabala S. S., et al., 2012, Monthly Notices of the Royal 
Astronomical Society, pp no-no 
Shapiro K. L., et al., 2010, Monthly Notices of the Royal 
Astronomical Society, 402, 2140 
Skelton R. E., Bell E. F., Somerville R. S., 2009, The As¬ 
trophysical Journal Letters, 699, L9 
Snyder G. F., Cox T. J., Hayward C. C., Hernquist L., 
Jonsson P., 2011, The Astrophysical Journal, 741, 77 
Springel V., Di Matteo T., Hernquist L., 2005, The Astro- 
physical Journal Letters, 620, L79 
Stoughton C., et al., 2002, The Astronomical Journal, 123, 
485 

Strateva I., et al., 2001, The Astronomical Journal, 122, 
1861 

Struck C., Smith B. J., 2012, Monthly Notices of the Royal 
Astronomical Society, 422, 2444 
Sturch L. K., Madore B. F., 2012, in TYPHOON observa¬ 
tions of the Lindsay-Shapley Ring, pp 180-182 
Thilker D. A., et al., 2010, The Astrophysical Journal Let- 


© 0000 RAS, MNRAS 000, 000-000 



20 T. Haines et al. 


ters, 714, L171 

Tojeiro R., et al., 2013, Monthly Notices of the Royal As¬ 
tronomical Society, 432, 359 

Toomre A., Toomre J., 1972, The Astrophysical Journal, 
178, 623 

Trager S. C., Faber S. M., Worthey G., Gonzalez J. J., 
2000a, The Astronomical Journal, 119, 1645 
Trager S. C., Faber S. M., Worthey G., Gonzalez J. J., 
2000b, The Astronomical Journal, 120, 165 
Trujillo I., et al., 2006, The Astrophysical Journal, 650, 18 
Trujillo I., Conselice C. J., Bundy K., Cooper M. C., Eisen- 
hardt P., Ellis R. S., 2007, Monthly Notices of the Royal 
Astronomical Society, 382, 109 
Vazdekis A., Sanchez-Blazquez P., Falcon-Barroso J., Cen- 
arro A. J., Beasley M. A., Cardiel N., Gorgas J., Peletier 

R. F., 2010, Monthly Notices of the Royal Astronomical 
Society, 404, 1639 

Weijmans A., et al., 2009, Monthly Notices of the Royal 
Astronomical Society, 398, 561 
White S. D. M., Rees M. J., 1978, Monthly Notices of the 
Royal Astronomical Society, 183, 341 
Wild V., Walcher C. J., Johansson P. H., Tresse L., Chariot 

S. , Polio A., Le Fevre O., de Ravel L., 2009, Monthly 
Notices of the Royal Astronomical Society, 395, 144 

Worthey G., Ottaviani D. L., 1997, The Astrophysical 
Journal Supplement Series, 111, 377 
Yang X., Mo H. J., van den Bosch F. C., Pasquali A., Li 
C., Barden M., 2007, The Astrophysical Journal, 671, 153 
Yi S. K., et al., 2005, The Astrophysical Journal Letters, 
619, Llll 

York D. G., et al., 2000, The Astronomical Journal, 120, 
1579 

van Dokkum P. G., 2001, Publications of the Astronomical 
Society of the Pacific, 113, 1420 
van Dokkum P. G., 2005, The Astronomical Journal, 130, 
2647 


APPENDIX A: CONSTRUCTING BURST 
MODELS 

A key prediction of the modern merger hypothesis is that 
merging galaxies with favourable gas dynamics undergo a 
centrally-concentrated burst of star formation. To determ¬ 
ine if the SFHs of our observed galaxies are consistent with 
this scenario, we construct a suite of synthetic SEDs that 
model the abrupt change in star formation during the mer¬ 
ger process. We begin with a continuous SFH model assum¬ 
ing an exponentially-decreasing star-formation rate (SFR), 
SFRbg ( t) = aoe~ t ^ T °, where ao is the SFR at t = 0 and to is 
the characteristic time. Atop this background of continuous 
star formation, we add Gaussians of fixed widths and amp¬ 
litudes to model the bursts during the merger process. To 
make our models as realistic as possible, we use two bursts: 
a small-amplitude, large width (i.e., long-lasting) burst rep¬ 
resenting star formation induced by tidal interaction, and a 
second burst with large-amplitude and small width repres¬ 
enting the starburst induced by coalescence (see Figure FATI) . 
The modern merger hypothesis also predicts the ignition of 
an active galactic nucleus (AGN) that provides a source of 
energetic feedback to quench star formation in the galaxy. 
To model this AGN, we impose an exponentially-decreasing 


SFR once the second burst has reached its peak to quench 
the star formation. 

Using ti for the time of tidal interaction and t c for the 
time of coalescence, the SFR due to both bursts is given by 

SFR b (f) = [Gi + G c ] g (t) 

with 


G x (t) = tti exp {— (t - t x ) 2 /2al} 

and 


9(t) 


1 t < t c 

exp {—t/n} t ^ t c 


where t x is one of ti or t c , a x is the burst amplitude at 
t = t x , (Tx is the width of the respective burst, and Tb is 
the characteristic time for the exponentially-declining SFR 
of the burst at coalescence. The fraction of stars formed (by 
mass) during the burst, Fburst, is defined as 


Fburst — 


stellar mass formed during burst 


total stellar mass 
Incorporating our definitions from above, this becomes 


Fburst — 


fto SFR b (t) dt 


f t f SFRb g (t) + SFR 6 (t) dt 


ft„ [G* + Gc] g (i) dt 


f f [a 0 exp {-t/r 0 } + Gi + G c ] g ( t ) dt 


(Al) 


where SFRbg is the star formation rate of the background 
model. We assume that this background star formation con¬ 
tinues during the burst such that the total number of stars 
formed is the background star formation plus the two bursts. 

Clearly, there are many parameters to be tuned in this 
formula. We choose to keep constant the widths of the 
bursts and the characteristic time of the background star 
formation. Specifically, we assume a total galaxy lifetime 
tf = t c + 3Gyrs, a characteristic time of tq = 2.5Gyrs for 
the background star formation, an interaction time of ti = 
t c — 250Myrs, an interaction burst width of ai = 750Myrs, 
and a coalescence burst width of a c = lOOMyrs. These para¬ 
meters are consisten t with the r esults of gas-rich major mer¬ 
ger simulations from lCox et al.l (l2008l f. Finally, we normalize 
the total (background plus bursts) star formation such that 
a mass of Mtot is formed. 

The amplitudes of the Gaussian bursts are independent 
of time, so we can rearrange the integrals in Equation [All 
using the substitutions 

rtf 

71 — / exp { - (t - U) 2 /2 gI } g (t) dt 

Jt 0 

rtf 

72 = / exp { - (t - t c ) 2 /2(7c } g (t) dt 

Jt Q 

/■*/ 

73 = / exp {—t/ro}g (t) dt 

Jt 0 


to generate a set of simple algebraic expressions to constrain 
the peak SFR parameters ao, ot Cl and ai such that 


a;7i + «o 72 = Mtot — 0073 

F burst 

ai 7 i + a c 72 = --—-0073 

J- -T burst 

If we assume the peak SFR of the interaction and coalescence 
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Figure Al. Star formation rates as a function of time (i.e., the 
star formation histories) for four of our constructed burst models 
used throughout the rest of this work: two in which 25% of the 
total galaxy mass is formed in a single burst (green) and two 
burst models with a 5% burst fraction (red). The solid lines refer 
to short e-folding times of 75 Myrs and the dashed lines to long 
e-folding times of 250 Myrs. The bursts are triggered at redshift 
0.12 when the galaxy is ~ 12 Gyrs old. Inset: zoom-in of star 
formation histories near the burst. 


bursts are tied such that on = ka c , then we arrive at our 
final expressions for the peak SFR parameters 

Mtot (1 Fburst) 

ao = - 

73 

_ Mtot Fburst 

Qc (£71 + 72)’ 

Using this model, we parametrize the star formation his¬ 
tory of our burst models using four parameters: the burst 
fraction Fb urs t, the relative strength of the peak SFR of the 
two bursts k, the characteristic time of the post-coalescence 
exponentially-decaying SFR. Tb, and the total mass formed. 
We assume k = 0.1 and Mtot = 10 11 Mq throughout, but 
generate our suite of burst models based on a variety of 
burst fractions and characteristic times. See Section [3] for 
these values. 


APPENDIX B: INDEX MODEL PLANE 

In Fig. [BU we present the detailed Balmer-D n (4000) index 
plane for each galaxy in our sample. For each galaxy, we 
show the suite of models outlined in Section [3] The inset in 
each figure shows the radial order of the index measurements 
with 1 indicating the central annulus. 


APPENDIX C: RADIAL SPECTRA 

In Fig. [CU we present the radially-binned, emission- 
subtracted spectra for each galaxy in our sample. The num¬ 
ber of spectra is determined by a signal-to-noise cut-off out¬ 
lined in Section [2.3.21 The spectra are listed top-to-bottom 
from the centre to the outer most annulus. For the central 
spectrum, we label several well-known spectral lines. Note, 
that due to the spectra being produced on different spectro¬ 
graphs with varying setups (see Section [2 .21 for details), the 
wavelength coverage is not the same for each galaxy. 
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Figure Bl. The stellar Balmer absorption versus 4000A break strength model plane showing two continuous SFH models (blue) and 
two burst models (green) covering a range of e-folding times permitted for each model type. The full suite of continuous SFH models are 
shown in light grey. Tick marks denote points along the model lines with known ages. For the continuous models, the ages are 1, 3, 5, 7, 
9, and 11 Gyrs. For the burst models, the ages are 0.2, 0.5, 0.9, 1.5, and 2.5 Gyrs. The four qualitative regions from the text are labelled 
and the large orange lines mark their approximate boundaries. At the top of each panel, we give the galaxy’s morphological classification 
next to its name. Inset: Zoom-in to the plane showing the indices marked in annulus order. 
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Figure Bl. (cont.) 
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Figure Cl. Observed spectrum (red) and emission-subtracted spectrum (grey) for each annulus with S/N/pixel > 10 from each galaxy 
in our sample. Each spectrum is scaled to arbitrary flux units to show detail and given in annulus order from top to bottom with the 
central annulus at the top. For the central annulus, several well-known spectral features such as the hydrogen Balmer series and the Ca 
II lines are highlighted with grey vertical lines and labelled accordingly. 
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Figure Cl. (cont.) 
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